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Generation of Microwave Power 
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P*-N-N™ structures have been fabricated in Si by boron diffusion into the n-type epitaxial layers supported by N* 
substrates. The current-voltage characteristics of these diodes have been measured beyond the avalanche breakdown |voltage. 
These characteristics have shown that these devices produce negative resistance. A source of microwave power in the X-band of 
frequencies has also been developed using these diodes in a tapered waveguide cavity. An average signal strength of ~ 12 dB 


has been obtained. 


1 Introdaction 

Since the historic discovery of transistor, there have 
been many attempts to develop semiconductor devices 
to operate in the microwave frequency range as 
amplifiers, oscillators, high spted switches, binary 
memories, etc. The ability of diodes and transistors to 
respond faithfully to high frequency external voltage 
variations is limited primarily by the minority carrier 
lifetime. Later, Au diffusion was used in semi- 
conductors to increase the recombination rate and to 
shorten the carrier lifertime by introducing deep 
energy levels in the forbidden energy gap. By this 
technique, the lifetime could not be reduced below the 


= nanosecond level and so p-n structures failed to 
a operate in the microwave range. 

= The observation of negative resistance region in the 
~~ forward current-voltage characteristic of tunnel diode 


: led to the possibilities of its applications in microwave 
S devices with such desirable features as: (1) negative 
me resistance at low dissipation levels, (2) the high 
frequency response, (3) low noise characteristics, and 
(4) relative invariance of its characteristics with 
ture. Being a majority carrier device, the 
tunnel diode has obvious advantages over transistors 
diodes due to short tunneling-time of carriers, but 
ts performance i is seriously restricted due to very low 
power handling capability and the difficulty of 
isolating the input from the output. 
Te iobeon generation of microwave frequencies using 
> biased p-n junction diodes aroused much 
tin the de characteristics of such diodes beyond 
4 gi psi voltage. Avalanche diodes 


structure is N*-P-J-P* and is called the Read diode 
after its discoverer. It consists of two regions: (1) a 
narrow avalanche region in the P-type layer in which 
the carriers multiply due to impact ionization when the 
device is reverse biased; and (2) a drift region in the I 
(intrinsic)-layer in which the carriers drift with 
saturated velocities due to high electric field. The 
negative resistance in the Read diode occurs due to the 
phase shift between terminal current and voltage of the 
diode. This phase shift consists of two parts: (1) phase 
delay of 90° caused by the avalanche multiplication: 
and (2) an additional phase delay caused by the finite 
transit time of avalanche generated carriers drifting 
through the I-layer. It is thus obvious that the external 
current is more than 90° out of phase with the voltage, 
causing the resistance to be negative. For a detailed 
study, the reader is referred to an article by Schroeder 
and Haddad’. 

The Read diode is a convenient structure for a 
theoretical analysis and physical understanding but at 
the same time more difficult to fabricate in the 
laboratory. The most widely used IMPATT diode 
structure is of the type P*-N-N* or its complementary 
N*-P-P* type. This three layer structure is’ much 


easier to fabricate than the four layer structure of Read 


diode. At high current densities, the current-voltage 
characteristics of these three layer structures reverse 
biased into avalanche breakdown region are 
determined by the space charge of current carriers. The 
influence of this charge can lead either to negative or 
positive static resistance depending upon the various 
device parameters”. The second mode of operation of 
avalanche diodes is called the TRAPATT mode—an 
acronym for TRAPped Avalanche and Transit Time’. 
For this mode of operation, the operating current 
densities are higher than for IMPATT mode but the 
frequency is much lower than the transit time mode 
frequency. In addition to these p-n junction devices, 
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there are some special materials from groups 3-5 
(GaAs, InP, etc.) which also show negative resistance 
but due to the intervalley electron transfer of electrons 
in the various conduction bands. For a detailed study 
of the utility of negative resistance in semiconductor, 
the reader is referred to the well known monograph by 
Roy’. | 

In the present paper, we have described the 
fabrication of various diodes of P*-N-N* type in Si 
and measured their current-voltage characteristics 
beyond the avalanche breakdown voltage. The dc 
negative resistance has been observed in_ these 
structures. A microwave cavity has been developed in 
the X-band using rectangular waveguide. These diodes 
have proved useful as good sources for microwave 
power when used in this cavity. 


2 Device Structure and Fabrication 

A typical diode structure and impurity profile are 
shown in Fig. 1 where V, = acceptor concentration, N, 
=donor concentration, NV, =surface concentration in 
the P* layer, N,=donor concentration in the 
substrate, N,=donor concentration in N-layer, x, 
= thickness of the n-type epitaxial layer and x,=P* 
diffusion depth. 

The fabrication process starts with a Si wafer 
consisting of N-type epitaxial layer supported by a 
heavily doped N* substrate. To form the Pt-N 
junction, boron was diffused into the epitaxial N-layer 
and the required junction depth x, and surface 
concentration NV, were obtained by subsequent ‘drive 
in’ process. This gives a complemeniary error function 
type of impurity distribution in the P* layer as shown 
in Fig. 1(b). Four-point probe method of resistivity 
measurement” was used to measure the surface 
concentration and impurity profile in the P* layer. 


i 
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Fig. I—(a) Device structure and (b) impurity profile 
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The junction depth and hence the width of the N- 
region after diffusion were determined by the method 
ion staining’. 

te diffusion, ch N* substrate was thinned down 
to remove any back diffusion and then high purity Al 
metal was used to make the ohmic contact to the 
substrate. The slice was then diced into individual 
diodes and each diode mounted on an Au plated Kovar 
pedestal with the N~ substrate in contact with the 
pedestal. A thin Au wire was then bonded on the P 
layer and then the standard Mesa structure was 
obtained to avoid the edge breakdown in the diodes. 
The diodes were closed with a Au plated Kovar cap, 
insulated from the pedestal by a ceramic ring housing 
the diode. 


3 Measurements 

To operate the diodes into avalanche region, high 
current densities are needed. If we want to go into 
IMPATT or TRAPATT mode of oscillations of the 
diode, the current densities may be much higher than 
this. Here we encounter one problem because if we 
apply pure dc bias to the diode, large currents will flow 
near the breakdown and the diode will be burnt due to 
excessive heating at the junction. To avoid this heating 
effect in diodes at high currents, the diode is operated 
under pulses of short duration superimposed on de just 
below the breakdown voltage. 

The experimental arrangement is shown in Fig. 2. 
Because of the short pulses, the current flows only fora 
short period of time and thus the excessive heating at 
the junction is avoided. The current through the diode 
/, is measured as the voltage drop (V, — V3) across a 
50 ohm resistance in series with the diode. The diode 
voltage V4 is just the voltage V, which was = 


directly on a calibrated oscilloscope. The pulse ay 


could be varied from 30ns to 1000ns and the 
maximum pulse amplitude was nearly 40 V. It was not 
Possible to take observations below 200 ns because of a 
little ringing problem in the output pulse. Alll the 
measurements were taken at a peak repetitior 
frequency of 1 kHz. The two capaci cuit 


tors in the circ 
serve the purpose of isolating the de signal from th 
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» Fig. 3—Details of the microwave cavity 


4 Microwave Cavity and Experimental Arrangement 
A Dolph-Chebyshev type of microwave cavity has 
been designed in a rectangular waveguide for operation 
in the X-band of frequencies (8.2 to 12.4 GHz). This 
type of filter gives the shortest length of taper for a 
given improvement in reflection coefficient, in 
addition to good impedance matching between the 
diode and the microwave circuit®. Such a cavity has 
been designed (Appendix A) and the details are shown 
in Fig. 3. The waveguide portion consists of a taper 
with a moveable short at one end as shown. To 
effectively dissipate the joule heat, an integrated heat 
sink has been provided below the diode. The bias to the 
diode is applied through the 50 ohm coaxial connector 
and a coaxial cavity which isolates the de and rf 


very 


> F 


ce > photograph of the complete experimental 
q * eocmaema for the oscillator is shown in Fig. 4. As 


Su. 


Bee: 


uae € 
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shown, the output from the cavity is fed to a waveguide 
type screw tuner and then to an attenuator. Finally the 
output from the attenuator is given to a crystal 
detector which is connected to a VSWR meter. The 
bias to the diode is applied as explained in Fig. 2. When 
a slotted section was inserted in the circuit (not shown 
in the photograph), the signal strength was measured 
on a VSWR meter as a function of the detector 
position on the slot. 


5 Results and Discussion 

The current (/,)-voltage (V4) characteristics of a 
diode are shown in Fig. 5 for various pulse widths w,. 
As the diode voltage is increased by increasing the 
pulse height, the diode current increases, reaches a 
peak at a critical voltage and then decreases. Thus 
negative resistance is produced at the device terminals, 
which is also found to be dependent upon the pulse 
width. If the de bias on the diode is increased, then the 
current through the diode also increases without 
changing the nature of the /,-V,, curves. After a certain 
voltage, the current /, again increases slowly with the 
voltage, giving a positive device resistance. 

In Fig. 6, we have plotted the strength of the 
microwave signal obtained from the oscillator circuit 
of Fig. 4, as a function of the detector position on the 
slotted section. The oscillator gives about 12 dB 
average signal strength, which should be sufficient for 
many applications. 

The same diodes could also be used for operation in 
the TRAPATT mode for which the operating current 
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chip bonding and the coaxial cavity designed for this 
oscillator but the experiments could not be completed 
with these structures due to nonavailability of time. 
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APPENDIX A 


Design of the Taper 

The relevant design equations for the cavity have 
been taken from Matthaei et al.’ If the waveguide is 
stepped from an impedance Zp to an impedance Zr, 
the impedance mismatch will produce a reflection 
coefficient 7, given by 


T, =41n(Z;/Z,), Zr>Zo =~ (1) 


For the Dolph-Chebyshev filter, the relation 
between the cut-off wavelength /,, and the reduction 
ratio P is 
P=T/T,=1/cosh(2nl/A,.) IF 
where / is the length of the taper. All wavelengths 
| shorter than /,, will give a reduction of mismatch equal * 
ea cucccs to or better than P. —- 

Let a and 6 represent the width and heigh 
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Fig. 5—Typical current-voltage characteristics of'a diode 
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Thus, we find that Z, =111 ohm and Z, =444 ohm. 
_ The minimum number of sections for the taper are 
ete to be 8, each having a length of 4.82 mm. The 


_ impednace along the taper is given by the relations. 

| W+x'N=(Z,/Z,)° /Z,Z3, sae (5) 
— M-x/)=(Zy/Z,)- CiZ,2.. 2>x> = < 

The value of Q has been given in Ref 7. Thus for 
various sections the impedance Z is calculated and, in 
turn, the waveguide dimension 6 is determined from 

_ Eq. (4). Calculated results are presented in Table |, 
_ where: x is the distance measured from the wide end of 
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Lead-lead oxide-lead tunnel junctions have been made in a cross-type geometry and the temperature dependence of the 
junction resistance has been studied in detail from room temperature to liquid helium temperature. High resistance yee 
(~ 1000Q) show an increase in their resistance by about 20 % as the temperature is lowered from room temperature (R ) . 
liquid helium temperature, while the low resistance junctions (< 1 Q) show an anomalous behaviour. The I- V characteristics ts 
such junctions indicated a negative slope at RT, which as the temperature is lowered increases, becomes infinite at a particular 
temperature and ultimately turns positive. It is shown that, even by adopting the four-probe technique for the measurement of 
conductivity, the electrode resistance makes a finite contribution and thus the effective resistance, instead of the true one, of the 


junction is actually measured. 


1 Introduction 
Tunnel junctions are very often made in a cross-type 
geometry by first depositing a metal film for base 
electrode, oxidizing its surface to form a suitable 
barrier and subsequently depositing another film 
perpendicular to the first one to serve as another 
electrode. The usual four-terminal technique as shown 
in Fig. 1 is used to measure the junction characteristics 
to elimmate resistances associated with contacts, 
measurement leads and the metal films serving as 
electrodes. It has been often observed! that the 
conductance of a metal-insulator-superconductor 
tunnel junction shifts discontinuously as _ the 
superconducting electrode turns normal. Osmun2 has 
recently shown that this discontinuous shift is due to 
the finite resistance of the electrode in the normal 
(nonsuperconducting) state which is not eliminated in 
the usual four-terminal configuration. While studying 
the temperature dependence? of the /- V characteristics 
of Pb-Pb,O,-Pb Josephson tunnel junctions whose 
electrode resistances in the normal state are usually 
comparable to that of the junction resistance, we have 
that it is almost impossible to measure the true 


V4 
tit 2 
~ tq 


jp MOTH, 


resistance of the junction in the normal state due to the 
effect of the finite electrode resistance. Only when the 
electrode film becomes superconducting, the true 
resistance of the junction could be measured. 


2 Experimental Details 
The Josephson junctions used in the present 
investigation were prepared on 1 in x lin microscope 
glass slides by thermal evaporation of 99.9999 °% pure 
lead in a vacuum of the order of 5x 107°Torr. The 
details of fabrication technique are described 
elsewhere* and only the salient features are given here. 
The films were made by evaporating the lead through 
metal masks. The junction dimensions were 0.1 x 0.2 
mm’. After depositing the first film, dry oxygen was 
admitted into the vacuum chamber through a needle 
valve and the substrates were heated by infrared — 
radiation to about 50°C by keeping an infrared lamp . 
Outside the glass belljar. The substrate rature | 5 
was monitored with the help of a copper ntan 
thermocouple. The substrates were, in general, 
maintained at elevated temperature for about lhrand 
were subsequently allowed to coolforthrinanoxygen _ 
atmosphere to fabricate low-resistance 
However, by suitably changing these parameter 
were able to make the tunnel junctions with 
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I-V characteristics of a low-resistance junction at 
different temperatures are shown in Fig. 2. At 4.2 K, the 
lead is in the superconducting state, the J-V 
characteristic resembles that of an ideal Josephson 
tunnel junction. As the biasing current is increased, 
zero voltage supercurrent due to the de Josephson 
effect is observed up to a critica! current value (3.4 mA) 
above which switching to superconducting energy gap 
value takes place, and the usual quasi-particle 
tunnelling curve is followed. In the reverse path, the 
junction follows the quasi-particle curve. Fig. 3 shows 
the variation of the critical current with external 
magnetic field which is produced by passing current in 
a solenoid. A typical diffraction pattern type behaviour 
is obtained. The zero-voltage supercurrent is found to 
be absent in high-resistance junctions and only the 
quasi-particle tunnelling curve is observed. These 
results indicate that the junction barrier does not 
contain any metallic shorts and only tunnelling current 
is present®. 

When the resistances of both high- and low- 
resistance junctions were measured as a function of 
temperature, entirely different types of behaviour 
were observed. For  high-resistance junctions 
( ~ 1000 ohms), the measured value of its resistance (R,,,) 
increased by about 20% in going from RT down to 
LHT. However, for low-resistance junctions, peculiar 
behaviour was observed. As shown in Fig. 2, negative 
Seat “wf V curves was observed at RT which amounts 
pxsutance. As the temperature was lowered, 
value of R,, recat ane'by an 
1 Ps ai >. It was also observed that, by just 
film thickness (which amounts to a 
in the electrode resistance) and keeping all 
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rs unchanged, the temperature at which 
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CURRENT IN SOLENOID (mA) 
Fig. 3— Variation of the critical current with applied magnetic field 
which is produced by passing current in a solenoid 


Table 1—Temperature Dependence of R,,, R, and R, 


(Values of R,,, R, and R, in ohms) 


Temp. R,, R, R, 
K 

290 —0.210 3.66 18.33 
175 —0.120 212 10.39 
110 — 0,060 1.50 741 
75 —0,030 1.09 5.38 
40 + 0,003 0.63 3.11 
17 +0,.030 0.29 1.41 

4.2 +0.45 0 0 

(R,=R,) 


junctions, the electrode resistance is negligible as 
compared to the junction resistance; thus R,, is not 
affected significantly. However, for the low-resistance 
junctions, as shown in Table 1 the RT value of the 
electrode resistances R,(top) and R,(bottom) is 
significantly greater than that of the junction 
resistance. The electrode resistance thus seriously 
affects R,, and makes it even negative. As eA 
temperature is lowered, the resistance of the r 
ae Se ae u 
junction (R,) remains more or less ¢ 
effect on R,, decreases with lov 

' ‘is explains q alitati re 
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the voltage V,,, is measured from the opposite edge. As temperatures. However, it may be noted that, i 
2 case, the resistances of both the electrodes — 


the current flows along the face of the oxide-it pele es 
comparable to that of the tunnel junction and 


undergoes a potential drop due to the finite electrode unction a1 
resistance. resulting in a voltage profile along the oxide resultant effect of the potential gradients at bot | the 
face. V,, is actually the minimum voltage across the electrodes will have to be considered. It is, therefo ii. 
desirable to derive a generalized expression for the fee 
resistance on similar lines. i. Peete 


oxide barrier. When the electrode is superconducting, 
the voltage is constant across the barrier surface and = 
equal to V,,. As a result of this, the measured 


conductance G in the normal state is greater than the . 
conductance in the superconducting state. Osmun? has The authors are grateful to the Director, N; 


derived an expression for the conductance for this case. Physical Laboratory, New Delhi for pr 
If Go is the true conductance and G(V,,)is the measured facilities and to Dr (Mrs) Madhu Prasad 


conductance, then 
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The de resistivity of Na,O-SiO,-CaO glasses containing Ce**, Ce** and mixtures of Ce?* +Cu* and Ce** +Cu2?* ions 
was studied over the temperature range 303-623K. The results were interpreted within the framework of the polaron 
hopping conduction mechanism of Mott and Davies. Experimental results indicated that Ce** forms stronger bonding than 
Ce** and that the lack of hopping sites in Ce** -glass is responsible for its higher resistivity value. Also some sort of electronic 
mutual interaction between Ce**, Cu* and Cu’* in the same glass matrix is considered to be of more impeding effect on the 
conduction process. The increased glass resistivity due to relatively low gamma radiation was ascribed to the impeding effect of 
the scattering of charge carriers and their decreasing mobility, while the increase due to tempering effect was ascribed to glass 
compaction. The relatively decreased resistivity at high gamma radiation dose was ascribed to the significant increase of the 


gamma-induced defects and charge carriers. 


1 Introduction 

Most of the rare earth (RE) metals (with the 
exception of Eu, Yb and one phase of Ce) are 
considered to have three valence electrons and this is 
the configuration assumed in all band studies done 
to-date. The atomic 5d and 6s states lie close to one 
another in energy and are widely separated from the 4f 
level. One expects the 5d and 6s electrons to form the 
conduction bands in metals and hence contribute to 

the conduction processes independently of the 4/ 

electrons. The contribution to the conduction band- 

width from the 5d electrons is expected to be somewhat 
less than that of the 6s electrons since the spatial 
overlap of the Sd functions is somewhat less than that 
of the 6s functions. Note that this is expected to beara 
strong similarity to the situation found in transition 

: metals (TM). 

Whereas the 5d and 6s electrons are far away from 
the nucleus, the 4f electrons are tightly bound to the 
atom and so they do not overlap the neighbouring 
atoms appreciably. As a result, the 4/ electrons will 
form a very narrow energy band in the metal. The d 
electrons are similarly, but less strongly, confined and 
they form energy bands of narrow to moderate widths, 

because, being mostly localized, they can tunnel 
through the potential barrier. Their quasi-band nature 
accounts for the fact that the d band electrons exhibit 
both localized and itinerant properties’. 

The electronic band structure of metallic Ce in its 

two allotropic foc forms (7 and 2) has been extensively 
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studied by Waber and Switendick?. Interest in this 
metal arises particularly from the #-y phase transition 
(which is thought to be related to a change in 
occupation of the 4/ level), and from the associated 
anomalous behaviour of the resistivity for the metal 
(by itself or as an impurity in other metals*). Burr and 
Ehara* measured the susceptibility at high tempera- 
tures (through the melting point) and found no 
evidence of a 4f to Sd electron transfer. Since no 
proposed configuration for cerium metal was generally 
accepted, five different atomic configurations 
4f?~ *5d*6s? (with x=0, 0.5, 1.0, 1.5 and 2.0) were 
determined from relativistic HFS wavefunction’. 
Waber and Switendick? found a great sensitivity of the 
position of the 4f band to the assumed atomic 
configuration. 

The importance of interband mixing (or hy- 
bridization) has come to some prominence as a 
consequence of the anomalous behaviour of cerium 
metal and its alloys®. This behaviour is thought to be 
associated® with the proximity of the cerium 4/ level to 
the Fermi level, £;. Alloys with cerium impurities are 
magnetic when the 4/ level is below Er, and 
nonmagnetic when it is above Ey (Ref. 7). The 
closeness of E4, to E; makes interband effects large in 
the case of Ce. 

The transitions of rare earth atoms have been 
extensively studied because these atoms, as additives in 
glass, have laser applications®, while those of 3d-series 
have been the cause of pigmentation in glass since 
historical times. However, even today Ce remains the 
least well understood rare earth metal, in either 
polycrystalline or glassy matrices. 

The present work concerns with the role of Ce (as a 
multivalence RE element) in the conduction 
mechanism of Na,O-SiO,-CaO-Ce,O glasses. Also 
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Table 1—Remarks on the Glass Samples with Variols 
Amounts of CeO, and CuO (in g) per 100g 
alency of Ited at 
g O, CuO Valency of Me 
— i ——— —— 1723 K under 
Ce Cu conditions 
G-I 0.35 3+ — Reducing* 
G-Il 0.35 - 4+ — Normal s 
G-II] 0.35 0.1 3+ 1+ Reducing 
G-IV 0.35 0.1 4+ -2+ Normal 


* 10g of tartaric acid was added per 100 g of the glasses. 


. - Z 2 + 
the effect of mixing Ce? * with Cu** and/or Cu?* on 


the conductivity of the same glass matrix is 
investigated. 
2 Experimental Details 
2.1 Sample Preparation 
Glasses of the base composition: 


72Si0,.18Na,0.10CaO(wt %) were prepared using 
A R chemical grade materials. Sodium and calcium 
oxides were added as their respective carbonates. 
Washed white sand was used to introduce SiO). 
Copper and cerium oxides were added as: CuO and 
CeO,. The batches of the glasses were melted at 
1723 K in Pt-Rh crucibles in an electrically heated 
furnace under the conditions illustrated in Table 1. 
The melts were cast into discs and annealed at 773K 
for 2 hr, after which the furnace was Switched off to 
cool to room temperature. The specimens were then 
ground and well polished to obtain discs 2 mm thick 
for electrical measurements. 


2.2 Irradiation Facilities and Electrical Measurements 

The glass specimens were exposed to different y-ray 
doses from a °°Co gamma source (6300 Ci) with a dose 
rate of 6.316 x 10* rad/hr at room temperature. The dc 
electrical resistivities of (the irradiated) specimens were 
measured after 2 hr of irradiation. A simple indirect 
method was used to measure the resistivity and the 
details of the electrical circuit are presented elsewhere’. 


The applied field used was 500 V/cm, which lies in the 
ohmic range. 


3 Results and Discussion 
Fig. | shows the dependence of In( p/T) wi 

p/T) with 1/T fo 

the _untempered and unirradiated glass pani 

Studied. This dependence is characterized by two 

temperature gions. The lower temperature region at 

which the resistivity increases with 


PURE & APPL PHYS, VOL 21, SEPTEMBER 1983 


region can be discussed in the following way. For a 
strongly disordered system (such as the glasses 
employed here), the main transport mechanism iS the 
electronic hopping, i.e. the change in the occupation of 
localized states, that are in contact witha phonon heat 
bath'®. Here a suitable formula for the resistivity, 
similar to that used by Austin and Mott"? is 


W 


RK 2 1 
in(o/T)=Ia| sree [+R . ( ) 
where R is the mean distance between the ions, C and 
(1—C) are the proportions of Ce** and Ce* 
respectively, v is the rate of decay of the wavefunction, 
(yy ~exp(— vR)) of an electron on Ce**, and 


w= Wit5Wo ee 

In Eq. (2), Wy is the energy of polaron formation 
and Wp the Miller-Abrahams term. Wp is also called 
the disorder energy, which is the energy difference 
between initial and final sites due to variations in the 
local arrangements of the RE ions. 

It is seen from Fig. 1 that glass G-II (containing 
Ce**) has a higher resistivity than that for glass G-I 
(containing Ce**). This can be interpreted if we 
consider that Ce** forms covalent bonds which are 
stronger than the Ce** ionic bonds. Thus it might be 
considered that the decrease in the hopping sites (and 
free charge carriers) in Ce** glass is responsible for its 
higher resistivity. 

According to Bandyopadhyay and Isard!2 it was 
believed that hopping between neighbouring ions of 
different TM elements is evidently more favourable, 
energetically, than the hopping between ions of the 


mae 


423 373 


623573 523 473 


323 303 
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same element at greater average distance. However, 
the results shown in Fig. | point to a different 
conclusion; the resistivity of G-III and G-IV is higher 
than that for G-I. This might be interpreted as follows. 
Ce** has 4 highly localized 4/-electrons. On the other 
hand, Cu** and Cu?” are of 3d itinerant electrons. 
According to Mott and Davis'* the hopping process is 
enhanced in the highly localized states, ie. RE- 
containing glass (sample G-I) must be of lower 
resistivity than of the RE-TM mixed one (samples G- 
II] and G-IV). This agrees with our present results. 
In a previous work®, the present authors 
investigated the conduction mechanism of sodium 
silicate glasses containing Cu, Cu* and/or Cu?*. 
Comparing these results with the present ones, it is 
found that glass containing Ce**, Cu* and Cu?* 
(sample G-IV) is of higher resistivity than glass samples 
containing Cu* and/or Cu?”. This may be attributed 
to some sort of electronic mutual interaction between 
RE and TM elements. This interaction is thought to be 
of impeding effect on the conduction mechanism. 
Fig. 2 shows the temperature dependence of In(p/T) 
for glasses from G-I to G-IV tempered at 873 K for 5 
hr. All tempered glasses are of higher resistivity than 
samples of the untempered ones. This result may be 
interpreted as that the untempered glass has an open 
structure of low density which permits larger fraction 
of charge carriers to participate in the hopping process. 
On the other hand, the tempered glass has the time to 
rearrange its structure to a compact one of higher 
density which inhibits the charge carrier motion. It is 
worth mentioning that gamma irradiation may also 


+ T,K 
623 573 523 473 423 373 


323 303 


temperature of In(p/T) for glasses 
tempered at 873K for 5 hr 


result in a glass compactness'*. This could be 
considered as one of the causes for the observed 
increase of resistivity by )-irradiation in the present 
work. The annealing effect, which has been observed in 
all glass samples, suggests that this effect is intrinsic to 
the bulk material and may be associated with bond 
rearrangement. This would support the assumption 
that the electrically active traps are associated with the 
bonding defects as suggested by Mott ef al.'° and 
Kastner e/ a/.'° 

The reciprocal temperature dependence of In(p/T) 
has been studied for all the glasses before and after 
irradiation with different doses of gamma radiation 
and similar results were obtained. Fig. 3 shows an 
example of this observed behaviour. From Fig. 3 and 
Table 2 it is clear that the value of the critical 
temperature 7, [defined here as the lower temperature 
limit of the In(p/7) versus (1/7) linear range] 
decreases with increasing radiation dose. This may 
suggest that the hopping process is enhanced at lower 
temperatures as the gamma radiation dose was 
increased as a result of the increase of the irradiation- 
induced defects. 

From Fig. 3 it can be noticed that the resistivity of 
the glass increases on increasing the gamma 
radiation dose. Fig. 4 shows this same effect for all the 
glasses studied, at 393K. From Fig. 4 it can be seen 
that glasses G-I and G-III are of similar beha- 
viour; the values of In (393 increase sharply up to a 
maximum value at about 0.19 Mrad. The resistivity of 
glasses G-II and G-IV decreases to a minimum at 
about 0.11 Mrad and increases to a maximum value at 


<T,K 


623 573 523 473 423 373 323 303 


16 1819 24 2-4 27 » OB 
(1000/T), K™ 
Fig. 3—The temperature dependence of In(p/T) for glass G-IV 
before and after irradiation with different doses of gamma radiation 
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Table 2— Val 


ey re os ee W (in eV) for All the Glass Samples Before and After Irradiation 
ues of /, 


Irradiated with a gamma 
Sample Unirrad- cas Gc MEER 
No. iated* 
é 0.09 0.13 0.19 0.25 
W ip W j W Fr. W f W 
0. 
378 0.86 371 0.90 367 0.96 360 0.94 352 93 
G-I : 
Hs 0.97 
G-Il a pe 359 0.81 355 0.82 353 0.87 328 
0.66 
G-III we va 363 0.85 348 0.89 346 0.88 340 0.95 
0.55 
G-IV ob oe 355 0.76 348 0.94 344 0.94 328 0.94 
(359) (0.65) 


* The values in the parenthesis are for untempered ones. 


Dose(Mrad ) 


Fig. 4—The effect of gamma radiation on the resistivity (p) at 393 K 
for all the glass samples studied 


0.3 


about 0.19 Mrad. Asa general observation (Fig. 4), at 
doses lower than 0.19 Mrad, the resistivity increases 
while for higher doses it shows a tendency to decreases. 

It is well known'? that most glasses possess 
submicroscopic crystalline structure in the form of 
helices or chains of silica. The substructure is basically 


Over distances up to 


correlation between adjacent SiO, tetrahedra. It is the 
substructure which forms the same electro 


bulk silica glass. Accordingly, the increase in resistivity 


B in their ility. 
addition, gamma irradiation has been Sahat 
cause compaction to the 
glasses subjected to pressures and was reco i 
resulting in a folding up of the SiO, tetrahedra rather 


than actual bond shortening'’. Such folding process, 
in addition to the rupture of bonds by gamma 
irradiation, '* is expected to increase the disorder of the 
glassy matrix. Consequently, this will result in an 
increase of Wp; the dwelling time also increases, giving 
a higher value of v. The dwelling time is defined’° as 
the mean time a carrier stays at a certain site including 
the possibility of leaving the site and returning. 

Thus, according to Eq. (1), the electrical resistivity 
should increase. The second and third terms of Eq. (1) 
can be proved to be the main factors for the increased 
resistivity of the y-irradiated glasses up to 0.19 Mrad. 
The value of the concentration of Ce?* ions (c) is 
expected to decrease according to the well known 
reaction!® 


Ce?* + hv>Ce*t +e 


If the relative concentration of the reduced state is 
always greater than that of the oxidized state of 
cerium, then the multiplication of (c) by (1 —c) will give 
products of larger values as the gamma radiation dose 
increases. This is because the amount of Ce** ions, i.e. 
(1 —c), will increase progressively. As a result, the value 
Rk/c(1—c) should decrease. Experimentally In(p/T) 
was found to increase (Fig. 4) up to 0.19 Mrad:; this 
indicates that the first term of Eq. (1) has small or no 
effect on increasing the resistivity in this rafige. It is 
reasonable that increasing the electron trappers (Ce**) 
will increase the impeding and scattering effects which 
will lead to the increase of the resistivity. 

At doses higher than 0.19 Mrad the number of the 
released charge carriers and defects caused by gamma 
irradiation increases significantly. Moreover, there is a 
Possibility for the reaction!®. aed 


Ce** +e Ce3+ + py 


to proceed in the highly- 
eoance the conduction and gh 


aN 


a 


B 
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Be -_sxlanton the decrease of resistivity of the highly 


glasses after reaching a maximum value 
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Results of ultrasonic velocity measurements in binary systems containing n-butanol as the common component pes 
used to evaluate several parameters. The other components are cyclohexane, acetone, cyclohexanone, cyclohexano and 
chloroform. The trend of variation of these parameters at constant temperature but increasing butanol mole-fraction is 
indicative of the nature of molecular association in the mixtures. 


1 Introduction 

Ultrasonic velocity and absorption measurements 
have been used by many workers in the study of 
molecular interactions in liquids’, polymers? and 
electrolytic solutions’. In this paper, we have 
evaluated a number of acoustic parameters for five 
binary mixtures containing cyclohexane, acetone, 
cyclohexanone, cyclohexanol and chloroform with n- 
butanol as a common component. The calculations are 
based on velocity (u) values obtained experimentally by 
Prakash and Srivastava* at 308.15K as a function of 
increasing mole-fraction of n-butanol, x(Bu). The 
significance of these parameters and their mathemati- 
cal expressions are given in literature’. 

The molar sound velocity (R = Vu’) is known to be 
independent of temperature, pressure and con- 
centration variations for unassociated liquids and for 
non-interacting molecules. Here these variations are 
shown in Fig. 1. The free-length (L,) or the average 
distance travelled by sound waves between two 
molecules is another parameter studied here. The ratio 
(L;/L,)° of the experimental value (L*) and the ideal 
value (L;) determines the relative association (RA) 
given in Table |. The free-volume V;=(V—b), where 
V is the molar volume and b is the van der Waal 
constant, the coefficient of viscosity 1(= Mu/K V2!) 
va a stag time t=7,) are given in Table 2. 

“ ee", - 
relutioer B, ia, is obtained from the 


Excess parameters® given by the diffe 
1g tence betwee 
experimental (exp) and ideal mixture (im) values, aad 


Fort and Moore® show that there is usually ; 
ally increase in 
excess volume and excess compressibility if strong 
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interactions between the unlike molecules are asbent; 
otherwise these are negative. 


2 Results and Discussion 

Variation of molar sound velocity (R) with x(Bu) 
shows that R is not constant so that the AB- 
interactions in these mixtures are expected to be 
associative. This is borne out by the rising trends of the 
relative association (RA), with increase in x(Bu). The 
fact that RA in n-butanol+cyclohexanol is less than 
one shows that in the pure cyclohexanol molecules, the 
molecular association is weaker compared to that in n- 
butanol. A similar statement could be made in respect 
of cyclohexanone. In n-butanol mixtures with 
cyclohexane, acetone and chloroform. however, RA 
variation is more pronounced. 

The Fort and Moore? criterion also favours the 
above conclusion. The excess adiabatic compressibility 
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Table |—Free-length and Relative Association as Functions of Mole-fraction (x) of n-Butanol at 308.15 K 


x 


u hx 10° E. A: RA u? (exp)/ 
m/s kg/m? m*/mol pm pm u* (im) 
| n-Butanol + cyclohexane 
0.000 1200 765.3 102.7 59.8 60.2 1.02 1.00 
0.540 1166 779.8 94.0 58.0 60.3 1.18 0.94 
0.700 1180 784.0 91.6 57.6 60.5 1.16 0.96 
1.000 1204 798.2 86.4 54.9 58.9 1.23 1.00 
II n-Butanol + acetone 
0.000 1123 773.6 69.1 60.7 63.7 1.16 1.00 
0.375 1133 790.9 74.9 56.1 62.6 1.38 0.93 
0.546 1164 798.6 77.5 55.1 60.7 1.42 0.99 
1.000 1204 798.1 86.4 54.9 58.9 1.23 1.00 
III n-Butanol + cyclohexanone 
0.000 1410 + 933.0 99.4 48.3 47.3 0.94 1.00 
0.403 1344 884.6 94.1 50.3 50.7 1.02 1.03 
0.628 1311 855.0 91.1 51.5 52:7 1.07 1.05 
0.771 1294 834.9 89.3 52.5 53.9 1.08 1.07 
[V_ -n-Butanol + cyclohexanol 
0.000 1440 937.9 101.1 49.1 56.3 0.84 1.00 
0.463 1374 880.0 94.3 50.5 49.8 0.93 1.07 
0.633 1329 852.8 92:1 53.0 52.1 0.94 1.06 
0.776 1298 53.8 1.07 


832.3 90.0 53.8 


V_ n-Butanol+chloroform 
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Table 2—Relaxation Time and % Deviations as Functions of Mole-fraction (x) of n-Butanol at K 
‘ Bs x 10"! V, x 10° n x 10* tx 10'5 (AV/V) Sri: ee 
m2/N m?/mol (N/m?) s s ifn 0 A 
| n-Butanol + cyclohexane 
0.000 90.74 7.19 0.633 57.4 0.00 0.00 0.00 
0.540 94.32 6.92 0.590 55.6 0.24 — 3.10 —0.78 
0.700 91.60 6.73 0.596 54.6 0.39 — 1.93 —0.24 | 
1.000 86.42 6.35 0.608 52:5 0.00 0.00 0.00 
II n-Butanol + acetone 
0.000 102.50 5.89 0.467 47.8 0.00 0.00 0.00 
0.375 98.49 6.11 0.508 49.9 — 0,87 — 1.80 — 1.58 
0.546 92.42 6.09 0.544 50.3 — 1.30 —0.28 — 1.50 
1.000 86.42 6.35 0.608 525 0.00 0.00 0.00 
III n-Butanol + cyclohexanone 
0.000 53.91 5.71 1.011 $4.5 0.00 0.00 0.00 
0.403 62.58 5.88 0.852 53.3 —0.22 127 0.11 
0.628 68.05 5.97 — 0.773 52.6 —0.31 2232 0.38 Ss 
0.771 71.53 6:01 0.727 52.0 —0.33 331 0.72 F 
IV an-Butanol+cyclohexanol 
0.000 51.42 5.65 1.061 54.6 0.00 0.00 0.00 
0.463 60.19 5.79 0.876 S27 0.23 3.15 0.72 
0.633 66.39 5.95 0.791 525 0.15 2.89 1.00 
0.776 Ales 6.03 0.731 S24 0.08 SAF 1.06 
V_n-Butanol + chloroform : 
0.000 75.71 5.51 0.836 63.3 0.00 0.00 0.00 | Sent : 
0.416 84.48 5.95 0.712 60.1 ~0.76 ~5.08 a 
0.586 87.26 gH 0.669 58.4 ~0.75 5.64 2.74 <n 
; : 304 -- 32148 >) ee 


0.740 8718  — 6.20 0.646 56.3 —0.52 ~4.04 1s: 


La 


the molecules’. In the other two systems, the u>(exp)/u?(im) shows similar ends 
associative effects are either very weak or absent. of n-butanol with | cy slo exane 
vay be Benes out here that the structuring chloroform are charact 
d here fe is not the same as that possible in interactions whi Ni 
€ hydratior of tert-butanol + water. This is of n-butanol mix! 
cial structure of water. However, a cyclohexanol, the 
tT Ir iia lather ) 
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cular interaction This indicates that the -OH group 

ni ol is unable to enhance the dipole moment of 

Pie cyclohexanone molecule. Perhaps their relative 

orientation is such that the resultant moment is 
_ decreased. 

The addition of a diluent to the polar and 

_ associating molecules of butanol may cause a 

_ progressive breaking of H-bonds of the latter if the 

petal is unable to form any complexes. This is not 

y to be the case for n-butanol + cyclohexanone 

n as the RA does not decrease. However, the two 

mm of atoms appear to interact very weakly. In 

1, the interaction between the n-butanol and 


ne molecules appears to be much stronger for the 
d above. The chloroform molecule has 


| | ies of this molecule ends to increase ihe 
free-length and adiabatic compre- 
dility. The bonds between the two types of molecules 


* 
>x perimenteé 


(in n-butanol+chloroform) are, therefore, not very 
strong. In the case of n-butanol + cyclohexanol, the 
positive excess compressibility and the increasing 
values of free-volume point to rather weak 
interactions. 
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The problem of calculating the quantum corrections to the thermodynamic properties of a mixture of molecules using hard- 


core potentials is studied. The explicit expressions are given for the first order correction to (i) the free-energy and pressure of a 
binary mixture of hard-sphere and square-well molecules and (ii) to the ‘excess’ free-energy of a binary mixture of square-well 


molecules. Numerical results are reported. The quantum effect, which increases with increase of density and with decrease of 


temperature is minimum at R= 1.0. 


1 Introduction 

This paper is concerned with the evaluation of 
thermodynamic properties of a semiclassical fluid 
mixture, whose molecules interact via hard-core plus 
attractive tail potential. We consider the case, when the 
constituent molecules have different diameters. For 
one-component fluid, this is usually calculated by the 
Wigner-Kirkwood (WK) method':?; when the 
potential is analytic and by the ‘modified’? WK 
method**, when the potential has a hard-core plus 
attractive-tail. This method can be extended to a 
mixture. 

In recent years, both experimental® and theoretical® 
attempts have been made to understand the structural 
and thermodynamic properties of a binary mixture of 
hard spheres and square-well molecules. However, 
these are confined to a classical fluid mixture, in which 
both species of molecules have the same diameter. The 
semiclassical fluid mixture has been studied less 
thoroughly. The only work available for a mixture of 
hard-sphere and Square-well molecules in the semi- 
classical limit is that of Singh and Sinha’:* who have 
considered the case when constituent molecules have 
same diameter. Very recently, Singh and Sinha®!° 
have Studied the hard-sphere mixture whose 
constituent molecules have different diameters in the 


mixtures whose molecules interact via a hard-core 
potential. The explicit expressions for the free-energy 
and equation of state are given there. We employ these 
expressions in Sec. 3 to obtain thermodynamic 
properties of a binary mixture (i) of hard-sphere and 
square-well molecules and (ii) of square-well mole- 
cular. The results are discussed in Sec. 4. 


2 General Theory for Semiclassical Fluid Mixture 
We consider an open system at volume V and 

temperature 7, made up of species 1, 2,..., v. There are 

N, molecules of species 1, N 2 Of species 2, etc., such that 


the total number of: molecules being N watt Ni. 


Further, we assume that the constituent molecules 
differ in size. In quantum statistical mechanics, the 
grand canonical partition function is given by 


¥ i 
- 


- ZN ‘ 
== Lier f- | Wall 2a, NY] | a 


mise 


we (2.1) 


zi=A; *exp[ Bui] . (2.2) 


where 4; is the chemical potential, A; = (2h? B/m)"/ is 
the thermal wavelength and. 


ge. Pe =| 1 aN 11,2, -+.,Nlexp . 4 


C—BAy]2.. ND (23) 
is the Slater sum of the System whose Se So 
given by: — 


Es 
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ujAl,m) = u? (I,m) + uP(I, m) .. (2.5) 


where u/(/,m) is the hard-sphere potential regarded as 
the reference potential and u/{/,m) is treated as 
perturbation, then the Slater sum can be expanded 
about the hard-sphere Slater sum* and the grand 
partition function can be written as’: 


mab wir | f-[mettanm 


x lite + fi? (I, m)] 
x mil ae + fii? (I, m, n)] ii dr, .. (2.6) 
where 
; 5? (I,m) = fF(1, m) 
+[1+f2(1, m)] fi) (1, m) .. (2.7a) 
42. 
fet = 2g em 
aN 2 
Pea”) | .. (2.7b) 
. 12.p? 
Sih? (Lm) = a7 [V url, m).Vub(l, n) 
: +V,quP(l,m).Vquh(m, n) 
+ VU R(I,n).V,uk(m, n)] sx @2s7c) 
= exh .. (2.7d) 


he thermodyaam properties can be calculatéd if 
n function is known. Thus the free-energy of 
nis — a 


case when the quantum effect is small and attractive 

potential is weak. Then Eq. (2.9a) reduces to: 

A? = sho,» X{X; || ether ee | 
The Rita of the hard-sphere mixture in the 

semiclassical limit, valid to the first order quantum 

correction, is given by® 


BAns _BAi 
NN 


TN v 
amr de: xixjgi)(di)dzAj; a (2.10) 
where Aj, is the free-energy of the classical hard-sphere 
fluid mixture, d;; and 4;; are, respectively, the diameter 
and thermal wavelength of the molecules of species i 
and j. 

We use the van der Waals, one-fluid (vdW1) theory?” 
of mixture to calculate the thermodynamic properties 
of the classical system. This theory approximates the 
properties of a mixture by those of a ie pure fluid 
with the parameters: 


d3 =). x;x; dj; a (2.11) 
tJ 
Eo =d)? Y xix 64543; 3: (2.12) 
tJ 
and° 
Xo =d,” 2 xed (2.13) 


cee ee m 
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To evaluate the integral appearing in Eq. (2.9a), we 
have used for the mixture RDF, the values obtained 
from the Oth order approximation in the conformal 
theory of Mo et al.!* In this approximation, the 
mixture RDF g(?* is equated to that of pure fluid at 


reduced density pd3, i.e. 


gi (r/di;, pd3) © g0°°(r/dij, Pad) 13 Lea 
dj; 
Thus At=B o> X\X alge a} If? .. (2.20) 
where /{}) = 2n( pd3) EB gi's(r*, pda) 
1 
x uP*(r*)r*7dr* 221) 


r*=r/d;, and uf*(r*)=uP(r/d;;)/€;; 
I} is the pure fluid integral at reduced density po = pd}. 


Its value ffepends upon the nature of the perturbation 
potential u?;*(r*). 


3 Binary Mixture 

We consider a binary mixture euciinn of one 
component of molecules of diameter d,, and other 
component of molecules of diameter d,,. For the 
_ unlike interaction, d,, and ¢,, are given by!® _ 


d,.=(d,, +d») iesteel) 
€12=C12(€11812)'? «« (3.2) 


where ¢,, is an adjustable parameter. Further ? yI8 
given by” 


4u=[0ir +4ay2y? + (3.3) 
tie oc qu obs yeaah for a classical har sphere 


7 


Once free-energy is known, the equation of state for 
the quantum mechanical hard sphere mixture is given 


by: 


ol . (3.11) 


=3./2 re 
We now apply the theory developed in Sec. 2 to 
binary mixtures consisting of (1) a component of hard-— 
sphere molecules and another of square-well 
molecules, and of (2) both components of “square-well 


molecules. 


where Pf > AEA 


3.1 Binary Mixture of Hard-Sphere and Square-Well Molecules 

We consider a binary mixture of hard-sphere 
(species 1) and square-well molecules (species 2) and __ 
assume that the interaction between unlike molecules _ 
is given by the hard-sphere interaction. Thus, 


ui{r)= 0, r<d,, = 4 
= — 6, dj,<r<ydj; j Z et ss ee ; 3 
=0, r> dj; G ee o1) - 

where €,, =€,,=0, ¢,, #0 and y= 1.5. This de a 


interesting both because of its simplicity and because = 
of its extreme values for €,/€22 and 12/€22. a & ag 


For this system, 1=1{2)=0 and the first order 
perturbation ot Ass =e ; rh sie 


i = : ie 3 
ce ee 
where I 


Pen ear 
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with n=1; where P>=pd3. The values for the 
parameters are given in Table 1. Thus, the free-energy 
of a binary mixture in the semiclassical limit, valid to 
the first order correction, is given by: 


A BAj, A 
i B + Be €22 AP*+ (3)aar- .. (3.17) 
and the equation of state is 
a" BP hs Be,, PsP*+(A/d)P4* .. (3.18) 
where Pi?*= =o) . 3.19) 


3.2 Binary Mixture of Square-Well Molecules 

We now consider a binary mixture of square-well 
(SW) molecules and assume that the interaction 
between the unlike molecules is also given by the 


- square-well interaction. 


va 


 BAg AAs 
~ 


The ‘excess’ free-energy for a SW mixture correct to 
the first order correction is given by: 


A 
$ +,/ 2nxx297'3(dy aed) 


_7* ae ae .. (3.20) 
4 : ee do hese fe = 
abe sn ve hard-sphere, gi's(d,>) is given by ; 
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3 
where Ape =2x,x,(42) id ~weis2R 
0 


and Af*= / mx x2gi3(d12)d3, .. (3.28) 
4 Results and Discussion 
4.1 Binary Mixture of Hard-Sphere and Square-Well Molecules 
Eqs (3.17) and (3.18) have been used to calculate the 
free-energy (with respect to the ideal gas) and equation 
of state of the binary mixture of hard-sphere and 
square-well molecules in the semiclassical limit. 
The classical values are correct to the first order of 
Be,,=T*~ '. Theresults are reported for a range of p* at 


T*= 1-40 


X= xF 050 
A/d =O-20 
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cere 
Table 2—Values of Square-Well Force Parameters 


INDIA 


15O = 1-40 
4 a. Pair d e/k A* 
a “i alae ee Ne-Ne 2.382 19.5 0.9249 
oa A-A 3.162 69.4 0.2632 
Ne-A YR F- 36.79 0.3553 


ee eS 


Fig. 3— Values of BP/p for a mixture of hard-sphere and square-well 
molecules as a function of p* for x, =x,=0.5 at T*=1.40 


T” = 2.47 
X= %=0-5 


Med =O-2 —— 


Fig. 5—Values of BA,/N for a mixture of peri: molecules asa 
function of p* for x; =x,=0.5_ 


* 
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ee 


force parameters”° of this model for Ne and A repor 
in Table 2. The values of A* for pure systems 
obtained from the relation A* =h/d(me)'/?. To obtain 
d;> and €12 for Ne-A, we use Eqs (3.1) and ( 3.2). 
simplicity, €,, is taken to be a va 
z 


evaluated from Ey. (3.25). 
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The axial splitting parameter D in dicalcium barium propionate doped with Mn?" has been studied from 139 to 597 K, and 
the temperature dependence of D of phase-I of the compound has been fitted to conform to a quadratic form in the range 265- 
597K. The effect of the phase transition is to cause a slight deviation from the quadratic form above the transition 
temperature, over a range of about 50 K. The discontinuous change in the D-tensor at the cubic+* orthorhombic transition at 
267K is analyzed in terms of an abrupt slowing of the propionate group flipping (between its two equilibrium positions), — 
accompanied by a static distortion. The cubic phase is seen to be the time-average of two equilibrium configurations 


corresponding to the low temperature phase. 


1 Introduction 

Dicalcium strontium propionate (DSP), dicalcium 
lead propionate (DLP) and dicalcium barium 
propionate (DBP) belong to the family of metal 
double-propionates. At ambient pressure, DBP 
undergoes two phase transitions from the room 
temperature cubic phase-I to phase-I] at 267K 
(presumed orthorhombic) and to phase-III at around 

203 K (symmetry unknown). 
In a previous study’, we have analyzed the ESR 
_ Spectra of Mn** -doped DBP in the crystal planes (111) 
and (001) and in powder samples, as a function of 
> temperature. The zero-field axial splitting D.exhibits a 
ae discontinuous change at 7;_,=267 K and a change in 
___ Slope of D versus Tat 7; _,=203K corresponding to 
 first- and second-order transitions respectively. The I 
_ © II transition also shows a change in the orientation 


of the D-tensor axes by about +15K in both planes. 


= oe In the present study, we have extended the earlier 
: inar Its on the. variation of D with 


to cover a wide temperature range. This 


stoichiometric proportions, by slow evaporation at 
313K. Mn** was doped in the crystals by adding 
MnCO; to the solutions. ESR spectra were taken of 
powder samples. The spectra were recorded on a 
Varian-E-109 X-band spectrometer with 100 kHz 
modulation. The temperature was varied between 139 
and 597 K, using a cold gas flow method (stability and 
resolution +1 K). a = 


Soa 


3 Results 


The axial splitting D is extracted from the ESR ~ < 
spectra in the same way as in our previous wo i ‘- 


results of the D versus T plot were f 
quadratic form of Eq.(1). The use of thi 


' Pinas - 
-he ae 
ie a i. ws 
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value Dyin at a temperature Ti, gi 
In our earlier work,tin the temy 
474K, the results of D ver 
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Table 1—Quadratic Temperature Dependence of Axial Splitting (D) in DBP 


[Fitting of temperature dependence of D in different regions of phase-I with 7, = lower temperature limit, and 7, = upper temperature 
limit. ] 


Case Region ¥; ,. Do 

(K) (K) (G) 
(a) @ 487 597 188.0 
(b) (ii) 374 474 235.5 
(c) (iii) 265 369 345.6 
(d) (i), (ii), (iii) 265 597 290.6 
(e) @, Gi 374 597 239.9 
(f) (ui), (iii) 265 474 324.3 


In the present study, up to 597K (Fig. 1), no such 
minimum is observed. We have re-analyzed our 
results, and with the new values of the parameters Do, « 
and B, Tmin lies well above the decomposition 
temperature of DBP (Table 1). 

We have analyzed the present data, assuming that 
the quadratic form of Eq.(1) will be closely followed 
far above 7; _), and that deviations from this form will 
occur near the transition temperature. If a single 
quadratic form is used for phase-I [case (d) in Table 1], 
the rms deviation A of experimental points from the 
: theoretical values [Eq.(1)] is seen to be large. 
Accordingly, we have divided the data for phase-I into 
three regions (i), (ii) and (iii) (Fig. 1) and least-square- 
; fitted the quadratic [Eq.(1)] to the data with the 
different combinations as shown in Table 1. 

The rms deviation A was seen to be different in the 
different regions, being largest in region (iii). 

In case (d), we also note that the value of Dpin 
exceeds the lowest experimentally observed value of D. 
In cases (c) and (d), the values of the parameters suggest 
minima in D, at high temperatures, that are not 
observed. Also, the deviation A,is larger than A,, while 
A, has the smallest value. Further, the coefficient of 
the quadratic term () increases in magnitude as the 
temperature approaches 7) 1 from above in sequence 
(i) (ii) (iii). The above observations support our 
contention that the effect of the phase transition at 

_ 7,_wis manifested in region (iii), and becomes larger as 


e the transition temperature is approached from above. 
One may question the use of three ‘regions’ in phase-l, 

and expect that it can even be two or four regions. This 
choice jis convenient in that one can compare the 
of the quadratic close to 7; -1 with that far 
_y, and also with an intermediate region while 


4 sufficient number of data points to make a 
re fit meaningful. The fact that such a 
‘the temperature range is not entirely 
t is the correct thing to do is borne 


A = 

may = Bal 

ee iy v 
Cee oe 


aerae-< 
a A 
J 
. ’ 
_ 7 


= P G , 
(x 10 (x 10 KD | D;. rms 
=, (K) (G) deviation 
A 

— 1.4355 0.67 10.69.4 43.9 0.0919 
— 1.8650 7 795.6 60.8 0.2677 
— 2.890 2.85 506.8 92.5 0.2832 
— 2.332 1.82 641.6 73.2 1.0726 
— 1.9208 1.26 762.8 64.1 0.2482 
— 2.6799 2.44 548.6 85.9 0.5109 


Eq.(1) is valid only far away from the transition 
[region (iii)] and as one gets closer to the transition, 
departures from the expected behaviour occur, largest 
being in region (i) which is the closest to the transition 
and less so in region (ii). 

The quadratic form [Eq.(1)] fitting regions (i) and 
(ii) together (case e) was extrapolated into region (iii), 
and the value of D so obtained was subtracted from the 
experimental D value in region (iii) (Fig. 2). The 
difference AD is plotted against T in the inset (Fig. 2), 
representing the effect of the phase transition. The 
dependence is approximately linear, and we attribute it 
to temperature-dependent modes corresponding to 
phase-II. 


4 Discussion 
4.1 Temperature Dependence of D 


There are two contributions to the temperature 
variation of D: one implicit—due to thermal 
expansion and the other explicit—due to harmonic 
lattice vibrations. Walsh? considered the explicit effect 


S 

Q 

Oo 180} 

Z 

F "6 yon 326 

— Temperature T( K) 

a o- Experimental 
nD x. Extrapolated 
: 


264 «2% 284 296 306 316 cy Tae) ee ase 7% 
- TEMPERATURE (K) 
Fig. 2—Departure (AD), close to the transition, of experimental 
values of D from expected values obtained by extrapolating the 
curve in regions (i) and (ii) —case(e)—into region (iii). The inset 
shows the approximately linear dependence of A D on Tclose to Ty-u 
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and showed it to result in a temperature variation: 


hy (3 
D=Do(1+8 coth34°7) ( ) 
where v is the frequency of the lattice vibration, and Do 
the static value of D. In the high temperature 
approximation (iv <k 7) this reduces to 


D=D,(1+ BT) 4 


A form similar to Eq.(3) was also obtained by 
Serway* and Pfister ef a/.> Marshall et al.° have 
incorporated the effect of anharmonic potential wells 
to obtain a form of the type: 

D=Do+xeothx4"+ Beoth? x0" ie 

The implicit temperature variation was removed? ~> 
by performing a variation of D with pressure 
experiment to obtain 1/D(0D/CP), with the thermal 
expansion, the compressibility and the equation of 
state of the material also being used. Leclerc and 
Manoogian’ also used the form: 


p= ]) eS + } 2 ee ional 


—p? hy 
D1 +aT+BT2)+D, coth (347) 


<>« (6) 


A possible reason for a quadratic dependence of the 
static contribution to D on temperature is obtained 
from the point-charge model, where the dependence of 
D on temperature is obtained from the point-charge 
model" where Doc(I/r"). Because of thermal 
expansion, r=ry(1+y7), when the expression is 
expanded binomially, it leads to the quadratic 
_ dependence. For the vibrational contribution, D may 
have a quadratic dependence, from Eq. (5) in the high 
temperature approximation. Thus, the form of Eq. (1) 
may be obtained from both static and vibrational 
contributions. Separation of the two can only be done 
by means of high pressure ESR experiments. A more 
general form than that of Walsh?, applicable at very 
low temperatures, was obtained by Shrivastava’, using 
the electron-phonon interaction, but that form does 
not concern us here. 

The pronounced nonlinearity of D(7) in phase-I 
observed, makes Eq. (4) inapplicable, wlhen wae 

resort to Eq. (1). It is our contention that the deviation 
from the quadratic form increases as TT, _y, and this 
_ deviation AD is plotted in Fig. 2. 


42 Discontinuous Change of D at 7, , 


According to the X-ray study of 
Glazer?. in phase-I the Pree Sethcen marten, 


Fig. 3—Low-frequency hopping in phase-I of the propionate group 
between two equivalent sites (from Ref. 2) 


splitting, D in phase-I is the time average, both in 
direction and magnitude, of two corresponding D 
configurations (D, and D,) in phase-II. The direction- 
averaging has been reported in our previous work!. 
These configurations correspond to our observed 
results of angular displacements by + 15° in the (111) 
plane, and + 16° in the (001) plane, of the D tensor axes 

in phase-II from those in phase-I. That is, if we 
associate with the D tensor, a frequency f.,. (~1/D), 

and consider a residence time in the equilibrium con- 
figurations a and b, in phase-I, 1/t>f,,¢ and in phase-II, 
1/t<f,,3 assuming that the time required to go from 
configuration a to configuration b is negligible in both 
phases. Thus, in phase-I, the direction of the projection 

of D (in either plane) is the same as the coincidence 
direction in phase-II, because of the averaging between 

the two configurations above 7,_ 1 (Fig. 4). 

In terms of the magnitude of the D tensor in the two | 

phases, we can obtain the line position Hj (in phase-I) 

by assuming that it occurs at the average of the line 
positions in phase-II, corresponding to the two __ 
configurations: pane 


H,= f+ Ai 


for various angles of 0. The averaging is sho ig F 

Fig. 4. — 

ys -I equals the value of D at the c denc 
rection in phase-II. The principal | oe Bee 
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| 
Fig. 4—Model for averaging of D corresponding to the motion of 
. the propionate groups above the transition. D, = 1/2 Dy(3 cos? 0 —1) 
is plotted as a function of 6, the orientation with respect to the 
; magnetic field in planes defined by the two D directions in the 
respective phases. D, is the time-averaged value in phase-I of D, and 

D, of phase-II 


from the direction cosines previously reported’. From 
the powder spectra, Dy is 165.2G, and according to 
-Eq.(8), D; should be 134.6G while it is actually 
_ 149.5G. This difference, we infer, arises because the 
4 ‘first order I1— II transition is accompanied not only by 
x a ‘reezing of the propionate group flipping motion, but 
1 ‘so by a static distortion of the oxygen octahedra. 

: “It may be noted that we have invoked a dynamic 
: m del to account (partially) for the discontinuity in D 
at Tis, and attributed the discrepancy to the static 
io = no the oxygen octahedra. On the basis of the 
aa we could just as well have explained 
ity entirely in terms of a static 
= Or ‘of the oxygen octahedron, since the 
ate group flipping may be already static on the 

scale. However, various 'H NMR second 


= 
, 
~~ 
& 
Root 


ssults!°"!! including the study of the second 
s a function both of temperature and 
BF perme in our laboratory!” "show 


ihe transition. Our ESR 
so ie Sine gree: 
ns corresp onding to 


cr ® Ss 


ae] 


phase-II, which is one of the possibilities considered by 
Stadnicka and Glazer? (in terms of the space group, 
Fd3m of phase-I is considered as a ‘dynamic’ mixture 
of left and right-handed tetragonal space groups 
P432,2 and P4,2,2, corresponding to phase-IJ). 

Such dynamic models, involving averaged tensors in 
the plane with higher symmetry, have been used by 
Bjorkstam'* and Schirmer and Muller!*. According to 
Bjorkstam, the averaged EFG tensor g of a deuteron in 
a ‘close’ and a ‘far’ configuration in an O...D—O bond 
in KD,PO, is 


gate =, (9) 

This statement is analogous to the one made above 
for the line positions arising from the D-tensor, which 
transforms similarly as the EFG tensor [Eq. (8)]. Also, 
it is clear that the transformation of the tensors D and q 
is not a simple vector projection (multiplying by cos 9), 
but is of the form: 


3cos? 0-1 
io cos ) 


5} . (10) 


Further at 7,_1, the discontinuity in D may be 
accounted for entirely by static distortion of the 
oxygen octahedron or by also incorporating the 
freezing of the flipping (on the ESR time scale) of the 
propionate groups (at the transition temperature). The 
latter model seems to be more consistent with the 
model of successive arrests of the motion of the 
propionates proposed by Aleksandrova et a/.'' and 
with the low frequency flipping of the propionate 
group at room temperature reported by Stadnicka and 
Glazer’. 

It may be noted that the observation of symmetry 
change at the I-II transition has been made only by 
ESR methods. The 7, versus 7 results of Sundaram et 
al.'* show no change at the first-order I-II transition 
nor do the IR experiments of Seetharaman et al.'® 
(although both groups detect changes at the second- 
order II-III transition). Probably, propionate group 
motions connected with the I-II transition do not fall 
within IR or NMR time scales. X-ray analysis of 
phase-II has not been reported so far. 
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The temperature dependence of etching rates of Makrofol polycarbonate plastic detectors has been studied. In the method 
employed, '$O-ions of energy 2.50 MeV/N (N signifies nucleon) are used as energetic ions for track formation in the detector. 
The experimental results indicate that the track etch rate is dependent on the energy loss rate of energetic ions; whereas the 
track registration sensitivity increases, the cone angle of '§O-ion tracks decreases at higher etching temperatures. The 
maximum etched track length agrees closely with the theoretical range. From the calibration curve, the critical energy loss is 
determined and is found to be 5.70 MeV mg’ ‘cm?, which is in agreement with that reported in literature. 


1 Introduction 
The production of tracks by energetic ions in solid 
state nuclear track detectors (SSNTDs) and the 
subsequent revealing of these tracks by chemical 
etching is a widely used technique for detection and 
identification of ions. One of the most valuable features 
of SSNTDs is the existence of a well-defined detection 
threshold, i.e. a critical value of energy loss rate, 
(dE/dx)... above which radiation damages produced on 
the SSNTDs by incident particles can be made into 
visible tracks by a selective chemical etching’. The 
critical value (dE/dx),,;, is characteristic of detector 
material. 
The track etch rate (V,) along the trajectory of a 
particle in the track detector is an increasing function 
of energy deposited around the track. For all cases of 
particle identification with nuclear track detector, it is 
necessary to know the relation between V, and the 
energy deposited along the trajectory of the particle 
under a certain etching condition. The curve of the 
etching rate as a function of the deposited energy is 
often called as the calibration curve. 
The track length provides the range of the particle 
and hence its energy. The maximum etched track 
length is also a measure of the mass A and charge Z of 
the incident particle”. The bulk etching at the detector 
surface shortens the lengths of the tracks. To correct 
the observed track length, it is essential to have an 
accurate knowledge of bulk etch rate (Vj). 
‘The present work was undertaken to draw the 
calibration curve of Makrofol-E polycarbonate plastic 
using 'O-ions of energy 2.50 MeV/N (N 
nucleon). The critical energy loss (dE/dx) ci for 
material is determined from the 
irve. The maximum etched track length is 
ith the theoretical range as well as with the 
by earlier workers. The dependence of 
emperature of etching is studied. 


+ 


— = 


2 Experimental Details 

Samples of Makrofol-E polycarbonate plastic 
detectors were irradiated with ‘§O-ions of energy 
2.50 MeV/N at an angle of 30° to the detector surface 
from cyclotron beams at the Joint Institute of Nuclear 
Research, Dubne, USSR. All etching was carried out in 
6N NaOH solution kept at a constant temperature in a 
thermostatic bath. The temperature at which etching 
took place was maintained constant to within +0.5°C. 
The temperature of the etchant varied from 50 to 90°C. 

The track lengths were measured with a transmitted 
light microscope ‘Olympus’ having an eyepiece 
micrometer of least count 0.215 um at a magnification 
of 900 x. The total error in measurements, arising 
from statistical errors, diffraction of light and 
microscope optical resolution, was found to be 
~0.5 pm. . 


3 Results and Discussion 
3.1 Effect of Temperature on Bulk Etch Rate V, and Track 
Etch Rate V, 

The bulk etch rate is measured by two techniques— 
(i) track diameter technique, and (ii) gravimetric 
technique. 

From the track diameter kinetics*** the relationship 
between track diameter (D) and V, can be expressed by 


p=2¥u pat | ae 


V+1 


where V=V,/V, and t is the etching time. 
If V>1, ie. V,>V,, then 


D=2V,t ow 


Eq. (2) is a linear relationship with a gradient of 2V,. 
Since V,/V,>1 for full energy fission fragments in 
plastics**, a plot of D of fission fragments versus 1 
should approximate to a straight line having a gradient 
of 2V,. 
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Samples of Makrofol-E are irradiated normally with 
a collimated beam of fission fragments from * Cf 
source in vacuum (~10°?Torr). The irradiated 
samples are then etched in 6N NaOH solution at a 
constant temperature. A graph is then plotted for 
Makrofol by taking the average diameter of about 250 
fission fragment tracks at various etching times. These 
plots are drawn for different etching temperatures. 
From the graphs, the values of V, are determined for 
different etching temperatures. 

In gravimetric technique*, the weight of the 
rectangular plastic sample (foil) before and after 
etching is measured with a microbalance. The bulk 


etch rate is given by 


_ Am 
~ 2Spt 


where Am (in g) is the dissolved mass in known etching 
time, t (in min) the etching time, S (in cm?) the surface 
area of the sample, and p (in g/cm”) the density of the 
plastic sample. 

The observed (projected) track length (/) is measured 
from the centre of the track ellipse at the etched surface 
to the end of the track tip. 

The corrected track length (L) is calculated by the 
relation? 

Vt 


rae ane Vilt—t,) 


V, 10*(in wm/min) < 19) 


. (4) 


where @ is the angle of incidence, V,t/sin 0 the surface 
etching correction, V,(t—t,) the overetching cor- 
rection and t, the time required to etch the tracks 
completely. 

The correlation between V, and L for an etching time 
t is given by the relation® 


L =| Ha 
0 


AL 
TN oe ihe | 


; Thus V, is measured by measuring the track length 
increase at different etching times. 


V, and V, are determined at 50, 60, 70. 80 x 
for 6N NaOH solutions. botenlad: 


The straight lines in Fig. 1 indicate that the 
dependence of V, and V, on etching temperature T 
follows the Arrhenius relationship of the forms 
V, = Aexp(—E,/kT) -- (6) 
V, = Bexp(—E,/kT) . (7) 
where k is the Boltzmann constant, A and B are 
constants, and F, and E, are the activation energy for 
bulk etching and track etching respectively. From the 
slopes of the straight lines, E, and E, have been 
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Fig. 2—L and | versus t (etching temperature, 70°C) 


calculated as E,=(0.72+0.08)eV and E, =(0.70 
+0.7) eV. The values are in good agreement with those 
reported by earlier workers’**. The ratio of E,/E, is 
found to be 0.97. ; 

The variations of | and L with t for the sample etched 
in 6N NaOH at 70°C are shown in Fig. 2. It can be seen 
from Fig. 2 that / starts decreasing after complete 
etching time (t,). This decrease is due to the fact that the 
bulk etching shortens the completely developed tracks. 
When bulk etching correction and overetching 
correction are made, L becomes constant beyond t.. 

The experimental values of V, at different L fot the 
sample etched in 6N NaOH at 70°C are shown in Fig. 
3. Similar plots can be drawn for 50, 60, 80 and 90°C. 
3.2 Dependence of Track Registration Sensitivity Vi= ¥,/V,) and 

Cone Angle (9) on Etching Temperature (T) niin n : * 
db relationship i wo 


Vt » Am/he 


FARID & SHARMA: CALIBRATION OF MACROFOL PLASTIC TRACK DETECTOR 
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Fig. 3— V, versus L (etching temperature, 70°C) 


(n(i/sin @) 


rounded terminal ends. The track lengths are 
calculated using Eq. (5). The track length distribution is 
shown in Fig. 5. The peak position of the Gaussian 
curve fitted to the historgram represents the most 
probable length. The most probable track length of 2.5 
MeV/N '8O-ion in Makrofol-E is <L) =45.50 ym. 

To calculate the theoretical range of '$O- ion in 
Makrofol-E, the stopping power and range equations 
of Mukherji and Nayak” have been used. The range (R, 
pai! of an ion in complex media is given by 

OE 

i) (dE/dx) Je OD 


where 


[(dE/dx)Je=4-YsLy.A{(E/dx)} 0] 


is the stopping power of the complex medium at the ion 
energy E (in MeV), {(dE/dx);}¢ the stopping power of 
the ith atomic species at the same ion energy, 


Aol =o y;A;) the molecular mass number of the 


medium, A,and Y; are the mass number and the number 
of ak a molecule Of the ith aes whieh the 
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Fig. 6—Variation of V, with the energy loss (dE/dx) 


stopping power remains constant, till the final ion 
energy Eo corresponding to the velocity Vo is 
reached’°. The range is then given by Eq. (10). Using 
the equations of Mukherji and Nayak, a computer 
program was made and with the help of the computer, 
the range and energy loss of '§O-ion in Makrofol- 
E(C,6H, 403), @=(1.20 g/cm*) were computed. The 
computer lists the energy loss and the penetration 
depth (i. range) for each energy until the energy 
reaches the value zero. The computed range of 2.50 
MeV/N '30-ion has been found to be R=48.28 ym. 
The maximum etched track length agrees closely 
(better than 6 %) with the theoretical range as well as 
with that reported by Tripier et al.'! 


3.4 Calibration Curve and Critical Energy Loss 

The (dE/dx) values corresponding to different track 
lengths (ie. penetration depths) have been obtained 
from computer output and a plot of these two has been 
drawn (not shown ). The variation of V, with track 
length is shown in Fig. 3. From these two figures, we 
deduce that V, is a function of energy loss rate (dE/dx), 
as shown in Fig. 6. The value of (dE/dx) at which V, 


es equal to V,1s taken as the critical energy loss 
ee a track etching below which no etchable 
track can be produced. The critical energy loss for 
Makrofol-E polycarbonate detector has been found to 
be 5.70 MeV mg” 'cm?. This value is in agreement with 


that reported in literature?"’’. 
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Properties of interacting bosons have been studied using the formalism of the theory of many-body systems. Results 
obtained for the isothermal compressibility show remarkable agreement with the experimental values on liquid *He. The mean 
square relative fluctuation in the number density of the system at temperatures below the lambda-transition temperature of 
liquid He* and the degree of correlation amongst the bosons at such temperatures are also studied using the standard methods 


of statistical mechanics. 


1 Introduction 
Calculations on the thermodynamic properties and 
the energy excitation spectrum of interacting bosons 
and comparison of such theoretical results with those 
of liquid *He had been a subject of a number of 
rigorous studies'~° in recent years. Of the various 
forms of the *He-*He interaction potential, the 
: Gaussian equivalent of the Lennard-Jones potential 
developed by Khanna and Das” could be quite 
successfully used in calculating the energy excitation 
spectrum for the interacting bosons. Kebukawa et al.* 
calculated the excitation energy by considering the 
phonon-phonon interaction to be through the Morse 
0tential. But their calculations led to a perceptible 
difference between theory and experiment in the 
extremely low momentum region and also in the higher 
momentum region above the roton dip in the 
excitation energy spectrum for liquid “He. 
Commenting on such discrepancies between theory 


eee 


d experiment Kebukawa et al. hinted* upon the 
essity of taking account of a long tail of attractive 

ial and also put forward the idea that the best 
> potential would be that of Lennard-Jones 
of the Morse potential. But the main difficulty 

ennard-Jones potential, in such a theory, was 
ability of its Fourier transform. Khanna 
wever, built the Gaussian equivalent of the 
es potential, having the required Fourier 
nd calculated the energy excitation 


ar 


a ~ 
be 


Using this potential from Ref. 2, we have, in this 
work, studied the connection between density 
fluctuations and spatial correlations of the bosons. We 


have used the expression for the space integral of the - 


correlation function C(r) in terms of the means square 
relative fluctuation in density. The isothermal 
compressibility 7; is then brought in to enable us to 
study the correlation function (or for that matter, the 
mean square relative fluctuation in density). The 
theory to obtain the isothermal compressibility for the 
system of bosons is then incorporated and an 
expression is derived for the isothermal compre- 
ssibility, involving, in the process, the ground state 
reaction matrix. We have seen that the mean square 
relative fluctuation can be written as kT times the 
isothermal compressibility, where k is the Boltzmann 
constant and 7 is the absolute temperature and thus we 
have studied the isothermal compressibility and the 
mean square relative fluctuation in density for the 
system of bosons interacting via a potential composed 
of a hard core plus two gaussians. 

The results obtained by us in our present 
calculations, for the isothermal compressibility bear 
very good agreement with the experimental values for 
liquid *He. In some of our earlier works (Ref. 6 for 
example), we have calculated the isothermal 
compressibility for the bosons assuming them to 
interact via a potential composed of a hard core 
followed by a square well. Comparing our present 
calculations with these earlier results we observe that 
better agreement between theory and experiment 
results, when the inter-particle potential is assumed to 


525 


re ic? alah vo ~? ee 
4 7 mI ‘ 
OO Pe ee — 


INDIAN J PURE & APPL PHYS, VOL 21, SEPTEMBER 1983 


be composed of a hard core followed by two gaussians, 
than when it is of the hard core followed by a square 
well type. Bet: 

As regards the mean square relative fluctuation in 
density we observe that in our present calculations, it 1s 
of the order of about one-tenth of 1/N, where N is the 
Avogadro number, at temperatures below about 1K; 
while it becomes of the order of 1/N at and above 7 
=1K. This result may be interpreted as indicating 
some sort of phase transition in the system around T 
=1K, corresponding thereby to the second phase 
transition observed experimentally in liquid *He at 
around 0.6K. 


2 Basic Theoretical Formulations 

The measure of the degree of spatial correlation in an 
interacting system of N particles enclosed in a volume 
V is given by the pair correction function C(r) defined 
by 


C(r)=n'(r)—1; 


where 77(r) is the pair distribution function. The 
product n*(r)d*r determines the probability of finding a 
particle in the volume element d*r around the point r 
when there is already a particle at the point r=0. The 
case of a classical gas composed of non-interacting 
particles is characterized by absence of spatial 
correlations and in such cases 7(r) is identically equal 
to unity and consequently the pair correlation function 
becomes identically equal to zero. 

Considering a macroscopic region V,, in the 
assembly, let there be N,, number of Particles in this 
region and defining the number N,,, which is fluctuating 
as a result of the switching on of the interaction by 


N 
Nm= ), Hr) Wes") 


F—fy = ra (1) 


where the function Lr) is such that H(r)=1 when r lies 
inside the region V,, and zero when r lies outside of it, 
ps mean square relative fluctuation in density is given 


<N2)—<N,>? ] 
PS Aa +0 j,,cowr] ak) 


where C(r) is the pair correlation function. 

Again, the space integral of the it correlati 
function C(r) is related to the mah small ina 
ssibility Z; through the well known relation (Ref. 6) 


| 
C(r)d?r kTy, ? «. (4) 
in density is given by 
a 
Pao [l+ okT yr 1] =kT yy i 
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To calculate 7, we proceed from the known relation 


ve tey 6 
Te PX ee 


To calculate 77 from Eq. (6) one has to know the 
pressure P. The equation of state is given by 


Po ee 

kT 1=27p \ r°dr 

x [exp{ — Vir)/kT} —1]g(r) ase a 

The radial distribution function (RDF), g(r), is 

related to the pair distribution function (PDF), 77(r), by 
the relation: 


2 
a ia os 


where p is the particle density. 
The PDF is the inverse Fourier transform of the 
structure factor S(k) and is given by 


n(r) = at | S(k) exp( — ik.r) d?k .. (8) 


This structure factor, in its turn, determines the 
excitation energy spectrum E(k) given by 
h2 k? j 
SK) = SFE) wes (9) 
This means that, when E(k) is known for an 
interacting system, Eq. (9) can be used to determine S(k) 
and then substituting for S(k) in Eq. (8), the PDF can be 


calculated to know the pressure P required for 


calculating x; from Eq. (6). This method of calculating 


the RDF from E(k) has been followed in our present 


calculations. 


Our earlier work! gives a method of calculating E(k) 


in terms of the matrix elements toooo Of the ground 
state reaction matrix, and we have: 


4m* 


h2k2 1/2 : ! 
E(k)= (sre) +7272 Notoo rd baie (10) 


which gives® ae . 


(2n)°h ar 


= 57) 2™* Notoooo| 42 mr/2 O<r< bear 
vs ate aia AS em 
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the ground state reaction matrix foo00 is given? by 
Ph? 
Notoooo = =7 8 g2 + 40° *p Volta — 13) Aes RS) 


where /ty and /1, are the repulsive and attractive ranges 
respectively after the hard core. 

Using Eqs (13), (11), (10) and (7) the equation of state 
for the presently considered system of bosons is given 
by 


P a 7 
eee be 2 plea — Hae) + ay/ MH H4 — Hr) 
_ 40? m*V (U4 — MpR)pa? 
3h 
82°/*m* Vo plu alte) 
h?RT 
x [(u dug) +ay/ m4 — Hrd) ... (14) 
From the above equation and using Eq. (6) one can 
calculate x7 from the following expression 


PET EV; 
5 — oar! a ees. 
x [(14 — MR) + ay/ M4 — Hr) 
82°'?m* pV oak Tus — HR) 
3h? 
162° ?m* pV6(H4 — HR) 
SS 


| -1 
x [M4 — MR) +ay/ nit, Hel} yi (45) 


When, however, the bosons are assumed to interact 
through a two-body potential composed of a hard core 
followed by a square well, i.e. a potential of the type, 


+o r<a 
un=}— Vo a<r<b vas (16) 
r>b 


the expression for 77 becomes 


1 A2kT 82m? pV ,a7kT 
a —. 


og, xm p VokTIb— a)°(b?— a2)! 
om — (17) 


rical Calculations and Results 
so (3) and (4), the values of x7 and F at 
atures T are evaluated, taking the 


A? = 33 

u2, = 0.2206A? 
uz =0.1103A? 
a =2.1A 


* 


m-=1.6 

m 

Vo=14.11 x 107! erg 
k =1.38x 107 '°erg/deg 


The reason for our choosing these values of the 
parameters lies in the fact that the potential considered 
in Our present calculations could very nicely bring in 
qualitative and quantitative agreement between the 
theoretical values of the excitation energy of the 
interacting bosons and the experimental energy 
excitation graph for liquid *He, as shown in Ref. 2. 
Table 1 gives the values of x7 and F (the mean square 
relative fluctuation in density) of a condensed system of 
interacting bosons for the potential composed of a 
hard core plus two gaussians at different temperatures 
below the 4-transition temperature of liquid *He. 

In Table 2 are recorded the values of x; versus T and 
F versus T, for a system of bosons interacting with a 
two-body potential composed of a hard core followed 
by a square well using Eqs (17), (3) and (4). In these 


Table 1—Variation of x, and F(=kT z,) with Temperature, 
for Hard Core + Two Gaussians-Type of Potential 


T, K ~710° ®cm?/dyn = F(=x7kT) x 1077? 


2.2 0.0108 0.0327 
2.0 0.0104 0.0287 
1.8 0.0100 0.0248 
1.6 0.0096 0.0213 
1.4 0.0093 0.0179 
1.2 0.0090 0.0149 
1.0 0.0087 0.0120 
0.8 0.0084 0.0002 
0.6 0.0082 0.0067 
0.4 0.0080 0.0044 
0.2 0.0077 0.0021 
ae ee See en WERE See TO 


Table 2—Variation of y7 and F(=k7 y,)-with Temperature 
for Hard Core + Square Well-Type of Potential 


T,K Xr, em?/dyn F(=kTy,z) 
22 —95.0x107'3 288.42 107 
20 —67.0x107'3 —184.92x10~!9 
18  —42.0x107'3 -—104.33x10"' 
16  —23.0%107'3  —50.78x 107-13 
14 —11.0%107'3 —21.25x10~'% 
12. —37.0x10-'* --61.27«10"'* 
10  ~80.0x107'S -—110.40x10~'S 
08 —74.0x10-'® -—81.70x10"'® 
06 —15.0x10"'? =—12.42x10-" 
0.4  —40.0x10724 —22.08x 10-74 
0.2 =—15,0x107°2. —414x10"% 
OS) py Se eee 
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calculations we use the following numerical values for 
the different parameters appearing in Eq. (17), 


Vo=14.11 x 107 '% erg; b=6A; a=2.1A; 
m*/m= 1.6; A? = 33 | 

These values were earlier used in our calculations op 
the thermodynamics of interacting boson system’’”’ 


4 Discussion 

Comparing our calculations shown in Fig. 1. with 
the experimental results? for liquid “He for its 
isothermal compressibility we see that very good 
qualitative agreement between theory and experiment 
is obtained in our present calculations. Further, from 
Table 1 it is seen that in our present calculations the 
mean square relative fluctuation in density for a 
condensed system of bosons increases as_ the 
temperature increases but nowhere in the temperature 
region 0.4K <7<2.2K, the density fluctuation 
becomes unity. This makes the system quite stable and 
also means that very few particles go to the excited 
levels from the ground state when a two-body 
potential, composed of a repulsive hard core followed 
by two gaussians is applied to the system. In other 
words, this predicts quite a strong correlation amongst 
the bosons in the ground state and also the presence of 
stronger attractive interactions amongst them. This 
result was also obtained earlier by Khanna et al.> when 
they studied the properties of interacting bosons using 
the same type of potential but with a different 
formalism. 

As regards the order of the fluctuation in density, we 
see that for the present potential the order of 
fluctuation in the number density is very much smaller 
than 1/N in the extremely low temperature region, 
where N is the Avogadro number, and it becomes of the 
order of 1/N only at and above about 1 K. This might 
mean the existence of a transition in phase of the 
assembly at about 1K, thereby corroborating the 
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0-4 0-8 1:2 1-6 2:0 2-4 
T(K) 
Fig. 1—Plot of isothermal compressibility 77 versus temperature (7) 


existence of a second phase transition in liquid He II at 
about 1 K, as was observed experimentally. 
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An attempt has been made to explain the behaviour of the effective thermal conductivity of two-phase gas-solid systems at 
different interstitial gas pressures, using kinetic theory. Dependence of gas and solid conduction upon porosity and particle size 
has also been estimated. Calculated values of the effective thermal conductivity have been compared with reported 
experimental results. On using the proposed theory in Russel’s expression and Hashin’s lower bound, we find a fair agreement 


with experimental results. 


1 Introduction 

Experimentally, variation of the effective thermal 
conductivity (K_) at different interstitial gas pressures 
of lunar materials'’* has been widely studied. The 
effective thermal conductivities of glass beads** and 
fibrous silica®~’ have recently been studied under 
changing interstitial gas pressure. In the references! ~’ 
cited above one finds that behaviour of interstitial gas 
has been treated qualitatively and some empirical 
expressions of Kg for silica®*® have been suggested. 
However, the available expressions lack the generality 
and quantitative character. In the present work, 
expressions for gas conduction and Kr of atwo-phase 
gas-solid system under changing interstitial gas 
pressure are suggested. The behaviour of solid 
conduction at reduced pressure is also explained. 


2 Theory 

The thermal conductivity (K,,) of an unconfined gas 
having mean free path 7, density ?9, mean molecular 
velocity u, specific heat c, and specific heats’ ratio ) ata 
gas pressure Pp is estimated*’? by 


K,,= 49; — SJucv-oPo ma 


In a gas-solid two-phase system the interstitial gas is 
mainly confined inside pores and it behaves differently 
than an unconfined gas. When pore diameter d is 
comparable to the mean free path /, at the interstitial 
gas pressure, P, the collisions.of gas molecules with the 
absorbed molecules at the pore surface occur 
frequently. The total path length for N collisions is 
thus affected by the presence of a pore. The 
- consideration of total probability of intermolecular 
collisions for all values of molecular displacement 

gives an expression'’®'!° for effective mean free path 
aa ee. ; 


aie address: Teacher fellow, Govt College, Dholpur 


(A.). The effective mean free path 4, as derived earlier is 
given by 


ig 
d 


ed 

The interstitial gas of mean free path /, thus 
becomes equivalent to a free gas with a mean free path 
/,. Thus, analogous to Eq. (1), the thermal conductivity 
of interstitial gas at the interstitial pressure P may be 
written as 


l - 
K,= 4? —S)ettrpe re 


Here it has been assumed that ), u, c, and 
temperature of interstitial gas are constant at different 
interstitial gas pressures. 

Gas conduction—A significant variation in Kp occurs 
in moderate pressure region where mean free path A is 
comparable to pore diameter d. The thermal 
conductivity K, of the gas, using Eqs (1) and (3), is 


. (2) 


oro 
From the kinetic theory, for a free gas 
40, Pt. 5 ... (4b) 
A po Po 
As such Eq. (4a) can be written as 
AP 
K-K,( 75) eee (5) 


Now, consider a gas pressure P,, when the mean free 
path of free gas 4, is equal to pore diameter d. Using 
Eqs (4b) and (5) 


4AP=ApPo=AP. ees (6) 

We thus have 
KKB] Pg 
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P+ 53 \p. 
“J 


and as A.=d, using Eq. (6) we have 


A.) 5 4 
A+d)P. \P.+P 


Therefore, 


P 
K,= Ku aep) eee (8) 


Eq. (8) predicts the behaviour of gas conduction at 
interstitial gas pressure near to P,. The characteristic 
pressure P. may be defined from Eq. (8) as P- P, K, 
—+K,,, i.e. the pressure where thermal conductivity of 
a pore-confined gas is reduced to half of its unconfined 
free gas value (K,,). In addition, since A and dare equal 
at P= P., the molecular density will be uniform inside 
all the pores and on an average there will be at least one 
molecule per pore. The gas molecules accumulate 
inside the pore at 2 < d, while a few of the pores are 
vacant at 4 > d. Therefore, molecular density is 
optimized at P = P.. As the convective heat transfer is 
proportional to molecular density, the convection 
effects are also optimized at this pressure (P.). The 
radiative heat transfer is minimum at this pressure for 
this mode is much more predominant in a vacated 
pore, whence A> d. This indicates that thermal 
conductivity measured at the characteristic pressure P, 
would be solely an outcome of conduction process. 

High pressure region—When the interstitial gas 
pressure is much larger than P., the mean free path A 
becomes much smaller than d. Using Eq. (2), the 
thermal conductivity of gas approaches its maximum 
value (K,.), ie. as 2 < d, 


Ad 
and 
K,>K,, ... (9) 


Low pressure region—When interstitial 
pressure is 
much lower than P., the mean free path 4 becomes 
much larger than d. Using Eq. (2), we have, as 1 < d, 


ed 


surface will always contribute to conduction 
Therefore, K, cannot be made zero a allowed by Ea 


Characteristic pressure—According to Eq. (6), the 
characteristic pressure is a function of pore size only. 
As mentioned above in Eq. (6) 

] 
a see 
and the pore size dis proportional to D, the solid grain 
size®. Therefore, we may write 
] 
P aa} 

Thus it is indicated that characteristic pressure P, of 
an interstitial gas in a powder should be larger than P- 
of the same gas in beads of a solid material. 

Solid conduction—The solid conduction in an 
evacuated two-phase porous system is due to grain to 
grain contact. This can be accounted through area of 
contact of a solid grain?’''. The area of contact is 
proportional to the number of solid grains and area of 
a single grain. When grains under consideration are 
spheres of diameter D and ¢, is the volume fraction of 
solid phase present in a two-phase system, we then 
have the number of grains x(y,/D?).and the area of a 
single grain oc D?. Hence solid conduction K, is given 
by 
K,, =p mat (11) 
where f is a proportionality constant depending upon 
the nature of the solid material. 

Radiative transfer—The second important me- 
chanism of heat transfer in evacuated porous systems 
is radiative transfer of heat within the pore. This 
contribution should be Proportional to the pore area 
and the number of pores. The pore size d is 
Proportional® to grain size D. Hence if Wy be the 
volume fraction of gas ina two-phase system, we have, 


number of pores oc “ 
and pore area x D? 
Hence radiative contribution to heat transfer is 


K,,=09,D ... (12) 


| 
‘ 
‘ 
7 


— disa Proportionality constant 


thermal conductivity at very low p . ' ‘ee wee . 


Ke=K,+ K,, Pant me eH: 
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Effective thermal conductivity—tIn the literature we 
find a number of expressions'* ~?° for the prediction 
of effective thermal conductivity (K_) of a two-phase 
system. These expressions are meant for the system 
where gas is at a constant pressure, generally at normal 
pressure. These expressions for Kg may be used under 
varying interstitial gas pressure provided the thermal 
conductivities of gas and solid phases are inferred at 
those gas pressures. It is relevant to assume here that 
the thermal conductivity of solid phase (K) is invariant 
at different interstitial gas pressures. On substituting 
the appropriate values of K,,, P and P. in Eqs (8), (9) 
and (10), the thermal conductivity of the gas phase (K,) 
at different interstitial gas pressure may be evaluated. 
The variation in Kg occurs only due to the variation in 
K,. Therefore, the effective thermal conductivity Kr of 
the two-phase porous system may now be evaluated by 
substituting K,, @ and K, at related interstitial gas 
pressure. Here the porosity ¢ is assumed constant and 
independent of gas pressure. However, the value of Kr 
at very low pressure ( P < P,) may be evaluated by using 
Eq. (13). In general, K¢ is a function of K,,.K,, P, P-and 
@. Using functional form, we may write 


Ke=f(K,, K,, P, P., 9”) see (14) 


The behaviour of Ke near characteristic pressure P. 
may be estimated by differentiating Eq. (14) with 
respect to P and ¢ respectively. Thus. we have, 


ess 
(+)... =f (K,, ) P: and ~) 


c 


={"(K,, sae ~) as aE at P=P, 
or 
a) =constant (if @ is constant) (15) 
P /p=P; gt! 
and 


( =/'(K, Ky P) 
: P= P. 


(when K, is a linear function of ¢) 


| or 
; (&) =f (K, +K,,) 
a =constant (if P. is constant) --- (16) 


15) and (16) are equivalent expressions. Eq. (15) 
tes that variation in Kr with respect to interstitial 

s ure is uniform near the characteristic pressure 
ven two-phase system, provided the porosity 
constant. Similarly, Eq. (16) implies that 


oP - 
A | 
ae 


3 Comparison with Experimental Values 

We have tried to find the suitability of various 
expressions for the effective thermal conductivity of 
gas-solid two-phase system. Table | shows a 
comparison of the calculated [using Eq. (8)] and 
reported values of Ky of glass beads (size 50 ym) at two 
different interstitial gas pressures(10 ~* and 10° mm of 
mercury). It is seen that Hashin-Shtrikman’s lower 
bound’* and Russel’s expressions’® fairly agree with 
the observed value of Ky of the glass beads (size 50 um) 
system. The effective thermal conductivities (K;) of 
fibrous silica (fibre diameter 1.2 um) and glass beads 
(size SO and 150 um) at varying interstitial gas 
pressures have been calculated using these two 
expressions'*’'? only. Predictions of Kg from 
Lichtnecker’s expression'® are also satisfactory in case 
of fibrous silica (fibre diameter 1.2 um). Therefore, this 
expression has also been considered for further 
calculations. Figs I-3 show the calculated values of Ke 
along with the reported experimental values at 
different interstitial gas pressures. 

It is evident from Figs 1-3 Lichtnecker’s 
expression'® predicts Ke values nearly two times the 
experimental value of Kg for 50 and 150 um sizes of 
glass bead systems. However, the same expression 


Table 1—Experimental & Calculated Values of Ke (in 
Wm _'K ‘') from Various Expressions Using Eq. (8) at Two 
Values of Interstitial Air Pressure for Glass Beads-Air 


System 

[Porosity = 0.42; K, = 1.04; K, = 0.026; T = 300K; and size, 
50 wm] 

Expressions Ky x 10? K; x 10? 


at 10-*mm of H, at 1000mm of A, 


Repo- Calcu- Repo- Calcu- 
rted lated rted lated 


Wiener's expression”® 0.5 58 13 60 
(For parallel resistors) 
Wiener’s expression”° do _0.005 do 6.5 


(for series resistors) 


Lichtnecker’s 

expression'° do 0.75 do 22.9 
Landau’s expression’’ do 20 do 36 
Kumar & Chaudhary’s 

expression’ ® do 30 do 37 
Maxwell’s expression'® do 49.5 do 51.5 
Hashin-Shtrikman’s 

lower bound'* do 0.4 do 11.7 
Brailsford and 
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Hashin’s lower bound 

1 r’s expr. ae et 
[eo] a Litchtnecke ele ee 
[e=*] Russel's expr. -* 


4 


[a } Hashin’s lower bound 
[e=0]} Litchineckerss expr. 


[oa} Russel’s expr. 


12 - 3 Experimental /” ~; 
n 4 
7 Experimental x / 
'E 
t 
£8 ge: 
x) = 
=x w 
w ve 
x 4 1 
fe) fe) - 
o §6o 1 10 100 1000 1072 10"! 1 10 Lo-9) 1000 
P(mm Hg) P(mm Hg) 
Fig. 1—Variation in Kg of glass beads (size 150 um and porosity Fig. 3— Variation in Kg of fibrous silica (fibre diameter 1.2 um and 
0.38) with change in interstitial air pressure [The characteristic porosity 0.95) with interstitial air pressure [The characteristic 
pressure is 6.5 mm of mercury. ] pressure is 1.2 mm of mercury.] ae 


G4: 


gives Kg very close to experimental value of Ke for . 
highly porous fibrous silica (fibre diameter 1.2 um). 


It is also indicated that Hashin-Shtrikman’s lower _ 
bound'* predicts Ke values with a maximum error of 


10% for the considered systems. It is further noted that a 


Hashin’s lower bound 
Eee] Litchtneckerss expr. 
[=] Russel’s expr. 


predictions by Russel’s expression!®, using — 
proposed expression for K, are closest to the obser 
values of Ke. The maximum error from this e: 
is 8%. S =e 
The variation in observed value of Kj 
characteristic pressure is uniform in all 
— is systems as predicted by Eqs (15) and (16). Th 
6 3+ 1 wo. Solid conduction: and: makina 
P(e Hg) coefficients B and 5 have been « 
(having only a var , 
Table 2. These a 


Kex 10? (Wm! kK!) 


hese a 


values may be applied to systems made of the same 
materials having different porosities and grain sizes. 


4 Results and Discussion 
In a two-phase system comprising air and glass 
(silica) the ratio K,/K, is 10 ~? at normal gas pressure. 
We find that only a few of the expressions of Kr of a 
two-phase system are capable of predicting Kp at low 
K,/K, value. The ratio K,/K, further decreases with the 
reduction in interstitial gas pressure [also obvious 
through Eqs (8)10)]. This further reduces the number 
_ Of useful expressions for Kg for the present purpose. 
_ All those expressions which are linear sum of thermal 
z conductivities of phases would not lead to observed 
values of K, at largely reduced pressures, because an 
Ba. appreciable change in the volume fraction of solid 
affects very little the experimental value of Kr at very 
Seal * hea (Figs 1-3). On the other hand, 
expressions which employ the product K,K, account 
. Therefore, Russel’s expression’? and 
‘Sh irikman’s lower bound'* under the 
z Eqs (8), (9) and (10) are found to be the only 


e in characteristic pressure with a 
article size (Figs 1-3) indicates that a two- 
consisting of microfined powdered 


a a\sis 2 = 
Pets poe 


ewe 


4055, ere 
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Thermally stimulated current in naphthalene electrets prepared in 
presence of thermal excitation (7,=23, 55, 60, 65 and 70°C) and 
magnetic field (H , = 14.50 kOe) is analyzed using three methods, viz. 
(a) Garlick and Gibson initial rise method, (b) Bucci plot method 
and (c) Cowell and Woods curve fitting method. Trap depths from 
1.024 to 1.206 eV and reciprocal of trap escape frequency (to) from 
1.52 x 10° '’ to 2.24 x 10° '* sec are obtained. The spread in the value 
of trap depth shows that the traps are distributed and closely spaced. 


Persistent internal polarization in naphthalene 
thermoelectret has been studied by Shrivastava and 
Ranjit Singh’. According to them the phenomenon is 
due to trapping of both holes and electrons. The decay 
of charge is supposed to be due to retrapping in shallow 
impurity levels. Campos et al.? have suggested the 
presence of pseudo dipoles to explain the thermo- 
electret effect in naphthalene. Thermally stimulated 
current (TSC) in naphthalene has been reported by 
Campos and Mergulhao’, and they have obtained the 
value of 1.26eV for activation energy for a TSC 
maximum at 38°C. They concluded that trapping is 
taking place in deep as well as in shallow traps. 

Surface charge decay of naphthalene magneto- 
electrets has been studied by Agrawal and Bhatnagar*. 
It is suggested that, trapping of charge carrier in 
structural defects and grain boundaries, and pseudo 
dipoles caused by nearby trapping of an electron-hole 
pair, may be the cause of the observed magnetoelectret 
State in naphthalene. 

Thermally stimulated current technique, being a 
powerful tool to study the internal mechanism, is used 
in the present communication for the first time to study 


naphthalene polarized in presence of magnetic fi 
and thermal excitation. pig 


Experimental Details 


contact. A series of magnetoelectrets were prepared at 
various temperatures (7, )(23, 55, 60, 65 and 70°C) and 
# magnetic field (H;) of 14.50 kOe by the usual 
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method*. Magnetic field was kept on for 1.5 hr, 1 hr of 
which was at constant forming temperature and half 
an hour during the cooling phase. 

Keithley electrometer 610 C was used to study the 
TSC thermogram. The arrangement for TSC was the 
same as used by Verma and Bhatnagar® and was 
enclosed in a chamber which was heated by Tempo 
Flexotherm heating tape. A linear heating rate of 
1°C/min was used to record the TSC. The whole 
experimental procedure was carried out in dark. Only 
red safe-light was used for taking observations. 

Thermally stimulated current is explained by the 
equation’ 


a E 
iT)=Aexo| ~~ Fes | exe( — ep AT | BS 
To 


where i(7) is the measured current, E the trap depth, k 
the Boltzmann constant, T the temperature in Kelvin, 
B the heating rate, t) the inverse of the trap escape 
frequency and A a constant which depends on trapping 
parameters. The experimental results can be analyzed 
using the following three methods: 
(a) Garlick and Gibson method® 
According to this method the low temperature tail of 
Eq. (1) is given by 
: E 
Inf)=C—7 st (2) 
where C is a constant. Slope of the Straight line plot of 
In (i) versus 1/T leads to the value of E/k and hence E. 
On differentiating Eg. (1) to get the maximum value 
of current, one obtains 
wit Kis 
0 BE exp(E/KT,,) + Q) 


The time constant for any other temperature can be 
obtained by the equation 
T=T) exp (E/kT) 


(b) Bucci plot method? 
In the derivation of Eq. (1) 


w. (4) 


traps at temperature 7. The numerator eeiaiipalar ss 
ig. 1. hown in 


NOTES 


ox 
- 
a 
= 
Ww 
- 


TIME (t),sec 


CURRENT [i),A 


TIME (t),Sec 
Fig. 1—Demonstration graph for Bucci plot method, i.e. plot of 
current i(t) versus time (t) [Inset time (t) versus temperature (T)] 


and Eq. (4) yields 


Int=Into + p= ... (6) 
Values of t at different temperatures can be obtained 
from Eq. (5). And slope of the Bucci plot between Int 
and 1/7 gives E and intercept on In t axis gives In Tp and 
hence To. 
(c) Cowell and Woods curve fitting method'® 
Substitution of x=E/kT and dx=—E/kT7dT in 
Eq. (1) and on subsequent approximation gives a 
simplified equation 


i(x) A exp[—x—Bexp(—x).x 7] ee AD 


When Eq. (7) is differentiated and equated to zero for 
the maximum in the curve, it gives 

_ (exp x*)(x*)° 
==, +2 

Here x* = E/kT,,,, T, = temperature at the maximum 
current value and B’ an approximate value of B 
obtained by replacing the approximation of Eq. (7) by 
an equality. 

Eqs (7) and (8) are used for the curve fitting technique. 
The value of 7,, is obtained from the experimental 
curve. The method described in (a) and (b) are used to 
calculate an approximate value of E with T,,,; this leads 
to an approximate value of x*. This value on 
substituting in Eq. (8) yields B’. Value of i(x) is then 
calculated from Eq. (7). A is used as an adjustable 

3 to normalize the theoretical curve with the 
experimental one at the maximum, that is 


A=ilx)/expL —x*—B’ exp(—x*).x*~ Ai .. (9) 


-__ Anyerror in the initial value of E chosen is indicated 
by a failure to obtain a good fit between the 
experimental and theoretical curve. If the fit is not good 
sw value of E is chosen and the procedure is 


.. (8) 


repeated, until an excellent fit is obtained. The value of 
E corresponding to this fit is taken to be the correct 
value of trap depth. Value of to can be calculated from 
the value of B, obtained with the help of the equation B 
= 1/t9.E/BT. The value of t at any other temperature is 
given by Eq. (4). 


Results and Discussion 

A representative TSC curve is shown in Fig. 2. For 
this sample, the forming temperature (7,) was 60°C 
and forming magnetic field (H,) was 14.50 kOe. 

A single peak has been observed at 40.0+0.5°C. The 
peak currents for the various samples are given in 
Table 1, which also gives the value of trap depth E and 
T> Obtained by Garlick and Gibson’s method. 

Straight line plot between In(i) and 1/7 of Garlick 
and Gibson’s initial rise method for the sample 7, 
= 60°C and H,=14.50 kOe is shown in Fig. 3. Trap 
depth for the various samples varies from 0.764 to 0.977 
eV and Ty varies from 3.23 x 10° 1° to 9.41 x 107 ‘4sec. 

Bucci plot between In(t) and 1/7 for the same sample 
(T=60°C and H , = 14.50 kOe) is shown in Fig. 4. For 
different samples the values of E and ty, are given in 
Table 1. 


© EXPTL POINTS 
— THEOR. CURVE 


T¢ = 60°C 
H¢ =14-50 kOe 


CURRENT (i (t)],10 2s 


2%. 30 35 ath. ao BOyshs (a0 
TEMP (T),c 
Fig. 2—Demonstration of Cowell and Woods curve fitting method, 
i.e. plot of thermally stimulated current i(1) versus temperature (7 ). * 
Experimental points and —, theoretical curve for: (a) E= 1.079 eV 
(b), E= 1.081 eV and (c), E=1.083 eV] 


-30-0 
-30:2 
30-4 


-30-6 


-30°8 
-310 
31:2 


Oa 325 330 335 


(10/7 ),K™ 
Fig. 3— In(i) versus | 7 plot of Garlick and Gibson 


535 


INDIAN J PURE & APPL PHYS, VOL 21, SEPTEMBER 1983 


3-10 3-20 3-30 
(103/T),K 


Fig. 4— Int versus 1/7 curve by Bucci plot method 


The trapping parameters obtained by Cowell and 
Woods method are given in Table 2. This table also 
includes value of 1399 and Tma,, the value of t 
corresponding to room temperature and_ the 
temperature at the maximum of the TSC peak is 
obtained. The curve fitting technique of Cowell and 
Woods is demonstrated in Fig. 2. As seen from F ig. 2 
when E is less than the value of trap depth 
corresponding to the experimental curve, the curve at 
the initial stage lies above the experimental one and 
vice-versa. 

Garlick and Gibson method can give cnly a rough 
estimate of trap depth (E) because very few 
experimental points are available for the plot of In (i) 
versus 1/T. The slope of straight line plot so obtained is 
quite ambiguous. The values of trap depths calculated 
by this method deviate a lot from the results obtained 
by other methods as shown in Tables 1 and 2. 

The Bucci plot method gives better result, as the 
whole TSC curve is utilized to get the straight line plot 
between In tand 1/7. The values of E so obtained from 
Bucci method are used in Cowell and Woods curve 
fitting method as described earlier. The value of trap 
depth (E) corresponding to the best fit is taken to be the 
correct value (Fig. 2b). Trap depth calculated by this 
method varies from 1.024 to 1.206 eV and To from 1.52 
x 10°" to 2.24 x 19-14 sec. 

This distribution of trap depths suggests that instead 
of a single trap level, there are a number of closely 
Spaced traps. The traps may be due to inhomogeneities 
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Table 1— Values of Peak Current, E and To for the Samples 
Prepared at 14.50 kOe 


T, Imax Garlick & Gibson Bucci plot method 
method 
C. 105%4 
E, eV To, SeC E,eV To, S€C 

23 1.12 0.800 899x10~'! 1.077 2.50x oe 
55 0.90 0.960 1.81x107~'> 1.209 1.13~x 107’ 
60 0.90 0.844 1.51x107'! 1.080  1.88x10~15 
65 0.72 0.764 3.23x107'° 1025 1.92x107~'* 
70 1.10 0.977 9.41x10~'* 1048 7.72x10715 


Table 2—Values of E, t,T399 and Tmax for the Samples 
Prepared at 14.50 kOe by Cowell & Woods Method 


T; E, eV To, sec T300, Sec Tmax,» SEC 
23 1.078 2.26 x 107 1° 3.03 x 10° 4.72 x 10? 
55 1.206 $5275 *% 2.91 x 10% 4.19 x 10 
60 1.081 1.81 x 10715 2.66 x 103 4.68 x 107 
65 1.024 2.24 =x 10-14 2.54 x 10° 4.98 x 10? 
70 1.052 5.39 x 107 !5 2.61 x 10° 4.81 x 10? 


SR ee —™ 


at the grain boundaries of the polycrystalline 
naphthalene. 
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Effects of Electric and Magnetic Fields (A, +jA,)K*+(B, +jB,)K?+(C,+jC,)=0 _... (1) 
on Propagation Characteristics of tank 
Microwaves in InSb 
1=- 
P N GUPTA*, M RAM & SK TOLPADI pth w aa 
Electronics & Radio Physics Laboratory, Department of Physics, ax ci 
Banaras Hindu University, Varanasi 221 005 a ” | 52 
ae pat 9 as 
Received 29 June 1982; revised received 25 November 1982 B,= ve | 2+ & foe AZ {i 2} 
Characteristics of electromagnetic waves, propagating through a Pe w? a v2 9 2+A? 
semiconductor plasma in the presence of external electrostatic and 49h @ ve 1+ Ae 
magnetostatic fields, have been investigated. Dispersion relation is 
obtained for the case in which the applied magnetic field is in the (, _@*[,/0,\* 1 _, 
direction of propagation and the electric field in a direction A ee} LAr 


perpendicular to the magnetic field. From the dispersion relation, 
| expressions for attenuation and phase constants have been derived. _ w*/o,\[3+42 (a1\? 1 
Sample calculations have been made for the case of InSb under ~~? eo pie Go Be ee 
different conditions and the results obtained are discussed. It is 
suggested that such investigations are useful in obtaining with 
information about transport properties of semiconductors. 


The study of the propagation of electromagnetic waves 
> Sees a semiconductor plasma has gained ag=neyp,, static electronic conductivity A,=1oBo 
e because of its possible application to study 


< the transport property of semiconductors. Pepeone 
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was done by Bok and Noziers' B=o4y, | az 
nd Vural? Recently, Pic and Ligeon* have - n=2102°m™? rr  — 
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NOTES 


Table 1—Parameters Used in the Calculation 
S. No. n 
‘ m™? 
2 x10!" 
2 x10?° 
1.3 x 10?! 
1.3 x 107! 
1.3 x 10?! 
1.3 x 10?! 
1.3 x 10?! 
so e= 1.602 x 10° '°C; ¢,= 18.7; 
lo ae” 10° °H/m; &) =8.854 x 107 !? F/m 


-y in the range 1 to 15GHz while the 
: n constant («) decreases with the increase of 


gf? t is Seacsionliy independent of carrier 
a. and that it varies Hea with 
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magnetic field. It is observed that the variation is very 
small. 

Effects of drift velocity and magnetic field have been 
also considered and it is found that phase constant is 
nearly unaflected by the drift velocity. This result 
agrees with the experimental result obtained by Pic 
and Ligeon®. Variation in phase constant is negligible 
as the magnetic field varies from 0.1 to 0.4T. 


Thus, the study of the phase and attenuation 
constants in the microwave region in the presence of 
the external static electric and magnetic fields can be 
used to explore the properties of semiconductors like 
InSb. 
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Liquid metals show the existence of an excess ultrasonic 
absorption over the classical theoretical value. We suggest a two- 
state thermodynamical model for liquid metals which explains this 
excess ultrasonic absorption observed in eight liquid metals. 


The excess ultrasonic absorption observed in liquid 
metals by the absorption measurements!’? has been 
assumed to be caused by the existence of the structural 
bulk viscosity as in associated liquids. An ultrasonic 
wave passing through the medium causes periodic 
changes of pressure due to which the molecules of the 
liquids fall into the vacancies in the lattice during the 
compression phase and return to their original 
position during the rarefaction phase. This process 
gives rise to the structural relaxation which is 
responsible for the loss of ultrasonic energy. 

A linear relationship between In( n/n.) and AS is 
observed in liquid metals as in the case of associated 
liquids*, where 7, is the bulk viscosity, 7, is the shear 
viscosity and AS is the entropy of melting. Also (n,/n,) 
is fairly independent of temperature around the 
melting point. This suggests a strong Coulomb 
interaction between the ions and the electrons which 
may cause them to behave like associated liquids This 
guides us to propose a two-state thermodynamical 


model for the acoustical energy loss in liquid metals. 
The existence of a double structure has already been 
confirmed by diffraction experiments’. 

The thermodynamic states of the atoms in the liquid 
under thermodynamic equilibrium can be characte- 
rized by the difference in the structural arrangements 
corresponding to the difference in the free-energy 
states. A fraction of the total number of atoms will be 
in state I and the remaining in state II, at a particular 
temperature and pressure. The acoustic waves passing 
through the medium disturb the population densities 
of both these states. The equilibrium will be established 
after the lapse of a finite time called the relaxation time 
t. The states I and II above correspond to the 
condensed state and the liquid state, respectively. 

Employing the Stokes and Kirchhoff’s relations for 
classical absorption and Hall’s theory*’® the following 
relation is obtained for the absorption of ultrasonic 
energy due to structural causes: 


at ode 
(7). - Cj Bo @ 


where f, is the relaxational compressibility, c, is the 
ultrasonic velocity of longitudinal waves in the liquid 
and fo is the compressibility in the liquid state at the 
concerned temperature and pressure. For the 
calculation of («/f)str. using Eq. (1), it is necessary to 
estimate f, and t independently of the ultrasonic 
measurements. A relaxation treatment is applied to the 
configurational component of compression to derive 
the expressions for f, and T: 


B, = VAV/V)?/2RT{1 +cosh(AF/RT)} one ae 


where V represents the average molar volume, AF the 
difference in free-energy between states I and II, R the 


Table 1—The Experimental and Theoretical Values (from the Proposed Model) of the Ultrasonic Absorption (a/f?) in Eight 
Liquid Metals 


Metal T(K) (AV/V) ¢ 


af? 


aw 


- 


ly 


NOTES 


eae aeihinn ET the temperature and (AV/V) the 


ae change in molar volume. 


_ The relaxation time is evaluated using Eyring’s 


- theory of reaction rates’ and is given as: 


=(Vn/ RIV +exp(AF/RT)] . (3) 


ames ang change in volume (AV/V) is estimated 
g the concept of open packed structure in the 

d state and the close packed structure for the 

id metals concerned. The details for evaluation of 


oa for liquid metals have been discussed 


ere®. With the knowledge of the basic 
eters and the equations discussed above, the 
cal values of (a/f7),.,. due to structural cases are 

i. The values agree fairly well with the 


experimental results’? for the eight liquid metals 
studied by us (Table 1) 


One of the authors, B V S Murthy, would like to 
thank H H Sri Shivarathri Rajendra Swamiji and 
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The all valence NDDO-MO method, after correction by S- 

function, is used to calculate the bond angle and bond length for the 
equilibrium configuration of some radicals and molecules. The 
results so obtained are in good agreement with the experimental 
values. Unlike in CNDO or INDO method no rotational averaging 
of repulsion integrals is done at the expense of information about the 
directional properties of orbitals. Hence, this method should be able 
to handle the problem of electronic structure in molecules possessing 
one or more lone pairs and also in molecules in which d-orbitals play 
vital role in bonding. 
The NDDO approximation was subjected to 
considerable theoretical analysis, but very few 
numerical calculations were attempted. Roby 
formulated a theoretical NDDO-MO scheme which 
did not seriously jeopardize the simplicity of Pople’s 
scheme’. We have presently calculated equilibrium 
geometry of a number of molecules using Roby’s 
formulation with a view to ascertaining its capabilities 
and limitations and to prove the effects of several 
refinements which have been suggested. To provide a 
consistent interpretation for the results, a minimum 
basis set of Slater-type orbitals, expanded in terms of 
three Gaussians each by a variational fit?, was 
employed in the calculation. Clementi-Raimondi 
exponents” were used for the orbital exponents of 
atoms of the first row of periodic table. The exponent 
for hydrogen was taken as 1.2. The integrals occurring 
in such computations can be evaluated to a high degree 
of accuracy with considerable ease using GTO basis‘. 
In the scheme formulated, the atomic orbital basis set 
{ W} is transformed into Lowdin symmetry orthogona- 
lized (OAO) set {2}, using the transformation‘, 


= sh 
ames 1" (I) 


ce ee to be solved over the {A} basis 
F°C*=C'E (2) 
Pr =H" +G (3) 
where, H* and G? are respective! ; 
\ y the one- and two- 
electron Parts of Fock matrix F*, The MO-coefficients 
C* and F may be expressed in terms of their non- 


orthogonal counterparts as follows: 

F4=§-t/2 Fst 

Crasitc on a 
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where S is overlap matrix over the non-orthogonal 
basis set. The OAO basis does not pose any difficult 
problem because it is easily evaluated over the non- 
orthogonal basis. However, Roby and Sinanoglu* have 
effected a simplification by proving the following 
interesting theorem. 

The electron repulsion in integral supermatrix over 
the OAO basis may be equated to that over the non- 
orthogonal basis in which the NDDO approximation 
has been made, provided the Ruedenberg expansion°® is 
sufficiently valid. In mathematical form 

G*=Gyppo - (6) 
1.€. (u4V 5/0605) = (Us V 4/0 cP) OaBOcp . (7) 

The truncated Ruedenberg expansion for a 
minimum basis set of atomic orbital is obviously far 
from correct. However, the source of error may easily 
be pinpointed in Roby’s scheme, and hence the method 
may be subjected to systematic refinements. 

The scale factors’ are employed for Coloumb 
repulsion integrals in the present work. They are 
evaluated using the S-expansion technique, correct 
through second order in overlap®. For one-centre 


Coulomb repulsion integrals, the corrections are given 
by 
B 
‘Vil= VAs SY ALS Abra ae 
+(S?\e (Vis — VARY) -» (8) 


‘VaP = Val + SHS* RAC AS — 1784 


+ LuSarge— VE 


Cc 
+ 2, ValiS7MRS (VAP — VAS) 
+ (SV (Vie? — VES]  @) 


The core-valence separation was carried out closely 
along the lines suggested by Lykos and Parr®. The total 
wavefunction is written in the form 
YW = fl (yore yval) .. (10) 


Wher rand v"' are antisymmetrized functions for 
tne core and valence electrons Tespectively. .of is an 
ator which interchanges 
electrons between the core and valence shells. The core 
orbitals are taken to be the original unpolarized atomic 
<a 


NOTES 


where E core = |W H.W dt 
Eva = Iv edt 

In Eqs (11) & (12) 

"gS 4 iy Re n.) 


. (12) 
. (13) 


ne ni 
m1 live . (14) 
and Mal (net 1 ...., no +n,) 
Nc + ny n, +n, j 
= (b) +4 
u=net 1 core 2 ae ae i oe : (15) 
with 


Hee = KH + Y (20) — kau] . (16) 


In Eqs (12-16) #y is the kinetic energy plus the 
nuclear attraction terms, the remaining term 
representing the Coulomb and exchange potentials 
from the core electrons. 

This procedure of core-valence separation is justified 
only if all valence orbitals are kept orthogonal to all 
core orbitals. In the present set of calculations, the 2s 
orbital has been made orthogonal to the 1s core on the 
Same atom by the Schmidt procedure. The valence 
orbitals on any atom have not been orthogonalized to 
the core on another atom but the overlap is expected to 
be too small to have any contribution. The 
conventional and more convenient way of treating the 
core as a point charge collapse at the nucleus has also 
been studied. 

The theoretical calculation of the equilibrium 
geometry for a molecule involves systematically 
minimizing the total energy of the system with respect 
to all independent internal displacement coordinates of 
the molecule. The binding energy of a molecule is then 
the difference between the total energy of a molecule at 
equilibrium geometry and the sum of the atomic 
energies of the component atoms. Actually the position 
of the absolute minimum in the total energy of the 
system is specified by the equilibrium geometry”. 

A systematic study has been made of molecules 
containing the atoms H, C, N and F with NDDO 
theory with only one or two polyvalent atoms (C or N). 
If we denote polyvalent atoms by symbol A and other 
atoms (H or F) by X and Y, the class of molecules 
considered is AX,, AXY and AY, molecules. The 
following constraint is placed on the nuclear 
configuration: 

AX,, AXY, AY,:C2, symmetry. 
We preferred here the Schoenflies notation'® 
expressing point symmetry. | 

The calculations were performed by starting with an 
initial guess of the nuclear configuration and varying 
individual parameters, bond angle in turn, until a 
‘minimum in the total energy was found. In all these 


Table | 


Calculated & Experimenta! Geometries fo1 AH,, 
AHF & AF, Molecules 
RAH, A RAF, A hae 
Cale Expl Calc Exptl Calc Exptl 

Singlet states 

CH, 1.14 1.12 105.0 103.2 

CHF 1.13 1.12 1.29 1.31 104.1 101.8 

CH, 1.11 - 1.30 1.30 103.8 104.9 

OH, 1.02 0.88 — - 104.5 104.3 

OHF 1.06 = 1.19 - 107.0 _ 

OF, — — 1.27 1.4] 105.9 103.3 
Doublet states 

NH, 1.05 1.02 — — 104.8 103.3 

NHF 1.07 — 1.24 — 106.6 - 

NF, — — 1.26 1.37 103.0 104.2 
Triplet states 

CH, 1.08 1.04 — — 160.0 180.0 

CF, a sees oF 1.33 — 129.8 - 


cases, bond angles were varied initially in steps of 1° 
and bond length in steps of 0.1 A. After one complete 
cycle through all parameters, the step sizes were 
decreased by a factor of 10 for a second cycle. 

For the carbon compounds CH, CHF and CF, all 
calculated bond lengths agree well with the 
experimental values, and the valence angles for the 
singlet states are also in good agreement. The bond 
angle in triplet state of CH, is correctly predicted to be 
larger than that in the singlet; however, the NDDO 
optimized calculations still lead to a bent triplet form 
rather than the linear form suggested on the basis of 
spectroscopic evidence. 

For the oxygen compounds, we may note that the 
experimentally observed geometry of water is well 
reproduced. However, the theory incorrectly predicts 
the FOF angle in F,O to be larger than in water. The 
NDDO gives small values for O—F bond length. The 
calculated angle in OHF is probably also too large for 
similar reasons. All these results are presented in 
Table 1. 
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A fresh study of vibrational molecular constants of nitrogen 
dioxide molecule has been carried out following Wilson’s F-G matrix 
method using a general quadratic potential function and kinetic 
constants. The complete vibrational analysis of NO, is presented 
and the results are briefly discussed. 


Analysis of absorption spectrum of nitrogen dioxide in 
the visible region has been difficult for many decades 
because of its complexity’. Several authors attempted 
to elucidate the complexity of NO, spectrum. Some 
groups investigated small details of the whole 
spectrum, like the rotational analysis of different 
vibronic bonds. Other groups measured fine and 
hyperfine structures of specific vibronic levels. But no 
satisfactory analysis of the coarse vibrational structure 
of the NO, is available in the literature. Recently, 
Bolduan and Jod1' used the matrix isolation technique 
to investigate the vibrational structure of the NO, 
molecule in the ground state as well as in the excited 
electronic state. They reported the complete 
vibrational frequencies of this molecule in rare gas 
solids Ar, Kr and Xe. Using the ground state 
frequencies in the present investigation, the molecular 
constants, viz. potential constants, compliance 


constants, vibrational mean amplitudes, Coriolis 
coupling constants, centrifugal distortion constants, 
inertia defect and thermodynamic functions for 
nitrogen dioxide have been studied. The structural 
parameters used in the present investigation have been 
taken from Ref. 2. The symmetry coordinates and the 
other details of the procedure of the calculations are 
the same as those given in the earlier papers? °. 

The evaluated molecular constants are reported in 
Table 1. The interesting observation in this molecule is 
that the major potential constants as well as the 
interaction potential constants are positive. The N-O 
potential constant obtained in the present work is in 
the expected range. Further, the values obtained by the 
present method agree satisfactorily with the values 
reported by Arakawa and Nielsen? (f,;=10.927; f,, 
= 2.038; f,=1.125; fi,=0.930). The compliance 
constants are evaluated by the Decius method®. They 
naturally exhibit trends opposite to those relating to 
the potential constants. It may be noticed that the 
bond-bond interaction and bond-angle interaction 
compliance constants are negative in this molecule. 
Using the potential constants in Cyvin’s secular 
equation’, the symmetrized mean square amplitude 
elements and hence the mean amplitudes are obtained 
at 298.16 K. It may be seen that the mean amplitude for 
the bonded as well as the non-bonded distances 
obtained in the present work is in the characteristic 
range for N—O vibration. | 


Table 1—Molecular Constants 


Kinetic constants Force constants (mdyn/A) 


Compliance constants (mdyn/A)~! 


a ea 


(10° 7g) 
(X)NO, (Ar)INO, (Kr)NO, (Xe)NO, (Ar)NO, (Kr)NO, (Xe)NO, 
Ke 14598 ff, 10.6942 10.7026 10.6480 pn 0.1007 On 
1007 0.1006 ~—S0,1009 
ku 09150 fig 2.3102 2.3290 2.9589 tee 00186 —0.0188  —0.0183 
ks 0.5448 ff, 10928 1.0896 1.0856 ny 0.9759 0.9784 09824 
a tz 04691 fg 0.9410 0.9382 0.9347 Nz 0.0706 -0.0704 —00712 re 
Coriolis coupling oer Centrifugal distortion constants (MHz) Mean amplitudes (A) at 298.16K ‘ — 
. “aie (Ar)NOz (Kr)NO, (XeJNO; —_—(Ar)NO, (Kr)NO, (XeJNO. 
13 ~0. -t 318.5032 319.1 20. 0.0387 0.0387. oe. 
ane ee an 320.6801 1, 0.0387 0.0387 ae § 
Re ME MM er 
Sh3 0.9110 Coie 1.7707 wre 1.7817 wieneunniailihe pase a pea 
(0.87121) (1.843) 0 Mor ars: i 
ge... a re 0.2599 yet oth cae 
1 fa — UOT ‘ : 
“Values in parentheses are the observed values. 
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NOTES 


The Coriolis toupling constants which may be 
determined experimentally can, however, be evaluated 
from a reliable set of potential constants and the values 
can be used for a detailed interpretation of the 
vibrational spectrum. The Coriolis coupling constants 
are obtained according to the vector method of Meal 
and Polo®. As expected, the value of (°, is negative. 
Further, ¢-values obey the sum rule: 

Ci3+C33,=1 
The Coriolis coupling constants when compared with 
the values available in the literature? agree 
satisfactorily. 

The centrifugal distortion constants evaluated by 
Kivelson and Wilson’s'® method are presented in 

Table 1 and they are also in the expected range. The 
present values of centrifugal distortion constants are 
compared with the observed values'’ in Table 1. The 
close agreement between the present values and the 

Observed values points out the significance of the 

method of kinetic constants. 
The calculated values of the inertia defect of nitrogen 

- dioxide molecule using Coriolis coupling constants 

-__ and centrifugal distortion constants are also given in 

} Table 1, along with the experimental values*’!?. It is 

interesting to note that the calculated value agrees 

_ quite well with the experimental value indicating the 

| Sanieanc of the present method. 

Lastly, the thermodynamic functions of nitrogen 

dioxide are calculated for the ideal gas state at 1 atm 
ressure assuming a rigid rotor harmonic oscillator 

. The effects of spin contribution, anharmonicity 

on: rigidity are neglected. The contribution to the 


vibrational part of these thermodynamic functions has 
been calculated from the tables given by Johnson et 
al'*. The values thus obtained form the thermody- 
namic functions are available with the authors. These 
constants are useful in the interpretation of the 
experimental results on thermodynamic properties. 

The authors are thankful to Prof S Sathikh, Director 
of the Madras Institute of Technology, for his 
encouragement and the facilities given to carry out this 
research work. One of the authors (K GR) is thankful to 
the CSIR, New Delhi, for the award of a junior research 
fellowship which enable him to pursue this 
investigation. 
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Levels of **Si were populated using *’Al(p, 7) reaction at E, = 2.05 
MeV. Decay from the resonance at 13.6 MeV in ?°Si was studied. 
The lifetimes of the levels were measured using DSA technique. The 
lifetime of 13.56 MeV state as 27 + 5 fs is being reported for the first 
time. The measured lifetimes have been compared with the previous 
experimental values and are found to be in good agreement. Further, 
the present experimental lifetimes are compared with the theoretical 
predictions using PW, KLS and ASDI interactions. 


There has been a considerable difficulty in 
understanding the spectrum of ?8Si. The nuclei at the 
beginning of the s-d shell exhibit characteristic 
rotational spectra (0*, 2*, 4*); however, this structure 
breaks down towards the middle of the shell. 
Therefore, the structure of 7°Si has been studied in 
recent years by various authors! ? by means of 
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27A1(p,y) ?®Si reaction to obtain information about the 
energy, spin, parity, lifetimes, etc. of various excited 
states. The recent experimental information about ?°Si 
has been compiled by Endt and van der Leun'®. A 
13.56 MeV level in ?*°Si corresponding to the proton 
resonance at 2.046 MeV was investigated by Kennedy 
et al.'1 They assigned 2*, 3*, 4 as probable spins to 
this level. However, no measurement of the lifetime of 
this level was reported. Therefore, we decided to 
investigate the lifetime of this level in the present study. 
The lifetimes of other levels are also compared with the 
theoretically calculated lifetimes to see the effectiveness 
of various interactions employed in the description of 
this nucleus. 

Experimental details—The reaction *’Al (p,7) 7°Si 
was used to populate states in **Si at proton 
bombarding energy of 2.05 MeV corresponding to the 
resonance at 13.56 MeV''. The proton beam was 
obtained from the Chandigarh cyclotron. The target 
was SkeV thick for 2 MeV protons. The )-rays were 
detected by a 50 cc Ge(Li) detector having resolution of 


1050-1150 1250 


NOTES 


= 2keV at 1.33 MeV. A typical spectrum is shown in 
Fig. 1. The numbers marked on the peaks in Fig. | 
correspond to the transitions shown in Table | and the 
corresponding level scheme is shown in Fig.2. The 
detector system was calibrated for energy and 
efficiency using )-rays from '°*Eu, '**Eu and °°Co. 
The drift in the detecting system was checked 
periodically using °°Co source and observing the )- 
rays from '°F (p, x) '°O reaction along with the well 
known background peaks. Centroid positions of y- 
rays’ full energy photo peaks were measured with the 
help of first moment method. 

Analysis—The lifetimes were calculated using DSA 
technique from the positions of the centroid using the 
equation! 


AE 


Fe) = —_—_—__—_—_—__—_—_—_——_- 
_ (cos 0, — COs 05)B. Eo 


where F(t) is the attenuation coefficient, AE the energy 
difference at angles 0, and 45, Eoo the energy of y-ray 
and B the Z component of the velocity of recoiling 7°Si 
atoms. 

Since the cross-section for the capture reaction was 
low, the spectrum at each angle was recorded for more 
than 20hr to get sufficient statistics on the relevant 
peak. 

From 13.56MeV level, there were three major 
transitions as shown in Fig. 3. Lifetimes were measured 
for all the three transitions. Experimental F(t) was 
calculated using Doppler shift equations and was 
found to be almost same in all y-rays from this level. 


Table 1—Transitions and y-ray Energies for Various Peaks 
Corresponding to Numbers shown in Fig. | 


Peak Transition Energy 
No. From — To (keV) 
1 6276.4-4617.7 1658.7 
2 1779.10 1779.1 
3 8944.5 6888.8 2055.7 
4 4617.7-1779.1 2838.6 
5 13555 10666.2 2889.4 
6 9315.8-+6276.4 3039.4 
7 13555 10417.0 3138 
8 13555-+10375 3180 
9 10417 + 6838.2 3578.2 
,¢ 10 8543.2-+4617.7 3925.5 
11 10417-+6276.4 4140.6 
Se ee eg 13555+9315.8 4239.0 
re 13 8944.8-+4617.8 4328.8 
6276.4-41779.1 ; 


13555-+8944.9 


E (keV) T(fs) 


J(Ref 1) 8 
= 13556 


(2,3,4) 


27-26 


55t6 


680+*70 


Fig.2—Decay scheme of ?°Si observed in the present work 
corresponding to the resonance at E,=2.046 MeV 


Theoretical attenuation factor was calculated using 
a computer code based on the nuclear stopping 
approximation of Blaugrund'* and the stopping 
theory of Lindhard et al.'* The target thickness was 
taken into account explicitly in calculating the 
theoretical values of Fit). The target was divided into 
10 layers of equal thickness and F(t) was calculated for 
each layer and the average value of Ft) was calculated 
from the values of F{t). A theoretical F(t) versus t 
curve was constructed and the values of lifetimes were 
extracted from these values using the experimental 
values of F(t). 
Results and discussion—Experimental lifetimes of 
various levels are shown in Table 2. The lifetime of the 
level at 13.56 MeV being reported for the first time is 
the weighted average of the values of lifetimes for the 
three transitions from this level as discussed above. The 
present value of lifetime of this level supports the spin 
assignment (2,3,4)* to this level'’. 

28S) is a nucleus having doubly closed proton and 
neutron shell, hence several shell model calculations of 
positive and negative parity states have been 
performed. Recently, two different interactions have 
been used for the calculations. The first is Preedom- 
Wildenthal'* (PW), which is based on an effective d-s 
interaction with two-body matrix elements adjusted to 
fit the 22Ne and ??Na data. The other is Kahana-Lee- 
Scott!® (KLS) interaction which is renormalized 
interaction based on a realistic bare NN potential and 
among the several renormalized d-s interactions. 
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Be 2 
300 @=120° 
13.56 6.89 
200 13.56 —= 4.62 
13.56 —= 1.78 

100 
2 
= 9 20.5 kev s 27.5 kev 36.2 kev, 
© 200 @=0 

100 

ask 1675 1685 2235 2245 2255 2945 2955 2965 2975 


CHANNEL NUMBER 
Fig. 3—Shift in the centroid position for 13.566.89, 13.564.62 and 13.56—1.78 transitions at 120° and 0° 


Table 2— Various Levels with Values of Excitation Energy and Lifetime 


Level E. 


Y 
keV keV i eRcine 
Experimental 

Present Ref. 8 
1779.1+0.1 1779.1 680+ 70 — 
4617.7+0.3 2838.6 55+8 40+7 
6276.4+0.4 4497.3 890 +90 — 
6888.8 +0.5 5109.7 67+10 45+9 
8543.2 +0.7 3925.5 3147 <5 
8944.9+0.6 2055.7 89+10 100+8 
9315.8+0.7 3039.4 16+4 <5 
10417.6+1.8 4141.2 22+6 27+8 
13555.6+3.1 11776.5 27+5 = 


Lifetime (fs) 
Theoretical 
Ref. 10 Pw's ASDI’’ KLS'® 
interaction. interaction interaction 

700 + 20 727 790 817 
63+6 55 57 51 
1200 + 120 1230 780 9150 
50+7 28.6 27 937 
18+5 18.6 13.5 9.5 
85+15 <= wr ae 
10+5 aba —_ = 
26+5 — pan i 


I 


The shell model calculations by Van Eijkern!” for 
this nucleus have been performed using adjusted 
surface 6 interactions (ASDI) in Ids), 2351/2 and 1d3),> 
configuration space. The calculated results from ASDI 
wavefunctions have been obtained with bare nucleon 
MI single-particle matrix elements and an isoscaler 
effective charge e°=2.15e: the isovector effective 
charge has been kept equal to unity. 


Table 2 also shows the values of lifetimes obtained in 
the present experiment along with the theoretically 
calculated lifetimes using the above three different 
approaches, viz. PW, ASDI and KLS. It is evident from 
Table 2 that the theoretical lifetimes obtained using 
PW and ASDI interactions are almost of the same 
order and they agree well with the experimental values. 
However, K LS interaction results are not supported by 
our experiment. The theoretical calculations of the 


several high energy levels are still missing; the 
calculations of these levels will be useful to verify the 
above conclusion. 
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The fission track technique has been used to estimate uranium 
concentration in some toothpastes manufactured in India. Uranium 
content varies from 0.91 to 3.56 ppm in sélected varieties of 
toothpastes. 


The solid state nuclear track detector (SSNTD) 
technique, because of its simplicity and applicability, 
has attracted many users in diverse fields such as 
anthropology, archaeology, biology, medical sciences 
and industrial technology’. The heaviest naturally 
occurring trace element, uranium, plays an important 
role in these investigations. The technique has already 
been exploited by various workers? ~> for the trace 
determination of uranium in semiconductors, water, 
milk powders, human blood, plants, soils, Indian 
Cigarettes (tobacco), tea leaves, portland cement, 
detergents and soaps, coal, flyash, steel, etc. In our 
laboratory, we have applied this technique for the trace 
determination of uranium in some toothpastes 
manufactured by various agencies in India. The aim of 
this particular work is to estimate the level of uranium 
content in toothpastes for the studies of radiation 
health hazards. 

We followed the external detector method as 
suggested by Fleischer er a/.! and Fisher® for bulk 
determination of uranium in homogeneous solids. The 
procedure is well-established. Toothpastes were dried 
in an oven at 200°C for 24hr. The dried powder 
(50 mg) was homogeneously mixed with 100 mg of 
methyl cellulose used as a binding material. The 
mixture was pressed into a pellet by a specially 
designed hand-pressing machine. These pellets were 
then covered on both sides with lexan plastic discs of 
the same diameter as the pellet. One such pellet was 
made of the standard glass dosimeter. The pellets 
covered with lexan discs were enclosed in an 
aluminium capsule and were irradiated with a ther- 
mal neutron dose of 5 x {0!5 (nvt) from the 
CIRUS Reactor of Bhabha Atomic Research Centre 
Trombay. After irradiation, the lexan discs were 
etched in 6N NaOH at 70°C for 30 min and were 
scanned for fission track density. A blank lexan 
detector was irradiated along with the samples and 
checked for background tracks. It was found almost 
devoid of background tracks. 


Table 1—Uranium Estimation in Some Indian Toothpastes 


[U Content of glass dosimeter = 20 ppm; thermal neutron dose =5 


x 101° (nvt)] 

Trade Manufacturing Ucon- 
name agency tent 
of the (ppm) 
tooth- 
paste 
Signal Hindustan Lever Ltd, Bombay 0.91 
Closeup Hindustan Lever Ltd, Bombay 1.04 
Promise Balsara Hygiene Products, Bombay 1.48 
Cosmo Hindustan, Rimmer, Delhi 1.78 
Colgate Colgate Palmolive Ltd, Bombay 1.82 
Flash Flash Laboratories (Pvt) Ltd, Bombay 2.41 
Forhans Geoffrey Manners & Co Ltd, Bombay 2.71 
Binaca Ciba-Geigy of India Ltd, Bombay 2.76 
Neem _Calchemico, The Calcutta Chemical Co Ltd, 

Calcutta 2.76 
Vicco Vicco Laboratories, Bombay 3.56 
vajradanti 


The uranium concentration was calculated using a 
simple relation for external detector method!: 


“=H )D® 


where the subscripts x and s stand for unknown and 
standard respectively, U, the uranium concentration, 
T the fission track density, J the isotopic abundance 
ratio of **°U and 738U and R the range of fission 
fragments in mg/cm’. The correction factor (RJ/R,) 
approaches unity for most of the silicate materials and 
plastics’. 

Toothpastes consist of various agents such as 
polishing agents or abrasive materials, humectants, 
binders, sweetening agents, flavouring agents, 
surfactants, flourides which act as a germicide and 
water. Materials used under these subheadings may 
vary from one manufacturer to another depending 
upon their formulation. Minerals, inorganic and 
organic matters present in the toothpastes are 
generally contaminated with minor amounts of | 
uranium. 

The uranium content of various too 
calculated by using Eq. (1) is given in Table 1. 
uranium content has been found to vary from 
3.56 ppm. Signal toothpaste ! 
Hindustan Lever Ltd., Bo 
uranium content of 0 
Viccovajradanti 
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NOTES 


i psesence of minerals derived from the 
yurces. The uranium content in the toothpastes 
le low and is not harmful for society. 
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The specific gravities of copper and cobalt complexes have been 
calculated from their X-ray chemical shifts. The calculated values of 
specific gravities of cobalt sulphide (CoS), cobalt disulphide (CoS,), 
cobalt oxalate (CoC,O,), cobalt sulphate (CoSO,.7H,O), copper 
oxide (CuO) and copper sulphate (CuSO,4.5H,O) are in fair 
agreement with their observed values. The specific gravities -of 
piperidinium pentamethylene dithiocarbonatocobalt (II) (Co-PPD); 
5,5’-thiodisalicylatocobalt (II) (Co-TDSA); 1-phenyltetrazoline-5- 
thionocobalt (II) (Co-PhTT); dimethylaminoethanethiolatocobalt 
(II) (Co-DMAET); a-mercaptopropionatocobalt (II) (Co-TLA); 
potassium deacyanocobalt (II) K,[Co(CN),0]; tetra-amino- 
carbonatocobalt (III) nitrate [Co(NH3),(CO;)JNO;;_ tetra- 
aminodinitrocobalt (III) nitrate [Co(NH;),(NO,),]NO;;_ tri- 
aminotrinitrocobalt (III) [Co(NH3)3(NO.)3]; 5,5’-thiodisalicylat- 
ocopper (II) (Cu-TDSA); piperidinium pentamethylene 
dithiocarbonatocopper (II) (Cu-PPD); diacetylacetonatocopper (II) 
hydride (Cu-AcAcH); diacetylacetonatocopper (II) chloride (Cu- 
AcAcCl); diacetylacetonatocopper (II) bromide (Cu-AcAcBr); 
diacetylacetonatocopper (II) nitrite (Cu-AcAcNO,); and a- 
mercaptopropionatocopper (II) (Cu-TLA) have been calculated for 
the first time. 


X-ray absorption spectroscopy has been effectively 
used’~° to determine bond length, effective ionic 
charge and coordination number of solids. Chemical 
Shifts of the X-ray absorption edges have been studied 
by several workers in compounds*~ 7, complexes®~ !° 
and spinels’ "'"? and various possible physico-chemical 
qualitative explanations have been put forward from 
time to time. Recently, the authors’ group'*"!* have 
calculated the chemical shifts with the help of plasmon 
shifts for several compounds. In the present note, we 
report the specific gravities of some cobalt and copper 
complexes using the experimental values of chemical 
shifts of X-ray absorption edges. The calculations are 
based on plasma oscillations theory of solids. 

In the X-ray absorption Process, the inner electron is 
ejected and goes to the first unoccupied state above the 
Fermi level giving rise to an absorption peak. The 
ejected electron produces an interaction'S:!© of 
suddenly created or annihilated core holes with 
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conduction electrons, in such a way that electrons 
oscillate collectively. In the Bohm and Pines’’ theory 
of plasma oscillations, the formation of energy loss and 
structural dependence has been ascribed to such 
interaction. Couchois!®, Lewis et al.1° and Borovskii 
and Shmidt?° have also shown such a possibility in X- 
ray absorption spectra. In this case the mean electronic 
neutrality is formed from free crystal electrons and 
from the lattice ions. Plasma can be stimulated to 
collective vibrations and fluctuations of charge density 
by means of a falling electron. The primary electrons 
suffer characteristic energy losses which correspond to 
natural frequency of collective density vibrations of 
free electrons or valence electrons. So it is quite 
plausible to explore such vibrations for the calculation 
of the density of compounds and complexes. 

On lines similar to our earlier!*** calculations, the 
‘plasmon shift’ or chemical shift of X-ray absorption 
edge of a compound or a complex is given by 


b6E=(h@,)c—(hw@,)u .. (1) 
where (fiw,)-=plasmon energy of compounds or 
complexes, (hw,)y =plasmon energy of metal and 

ho, =h(4nne?/m)'/? . 42) 


In Eq. (2) n is the density of free electrons, e and m are 
the electron charge and mass respectively. Further, for 
the electron density we get the following equation: 


_ OL 3 
n= =e aa t + . oe (3) 4 
where L is the Loshmidt number, ¢ is the density, Zis _ 


the valency (effective number of valence electrons 
taking part in plasma oscillations) and A is the atomic 
weight (molecular weight W in case of molecule or 
compound). Putting the numerical values for plasmon 
energy in Eq. (2), we get the expression: — 
ho, =28.8(Zo/W)'! tne 
Eq. (4) is valid for free electron-like metals but to a_ 
fairly good approximation, it can also be us fc . 
semiconductors and insulators?!. The plasmon er 

for dielectric is given by !%?2 


on ener 
Yan 
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Table 1—Calculated Specific Gravities of the Complexes 


Complex Molecular formula (5 Bors ZX W, (OQ catc (Core 
eV (Ref. 24) 
. Co-Complexes _ 
CoS 7.4141.0 (Ref. 25) 7 91.01 5.4229 5.45 
CoS, 18.99+1.0 (Ref. 25) 13 123.03 4.6465 4.268 
CoC,0, 5.47+0.3 (Ref. 26) 17 146.96 2.8928 3,021 
CoSO,.7H,0 10.4 Ref. 27) 87 281.12 1.8159 1.94 
Co-PPD _—[Co(C, ,H2N,S,)Cl,.(H0)] 4.5 (Ref.27) 115 412.331 «(1.0642 — 
Co-TDSA  H,[Co(C,4H,0,S).4H,0] 7.6 (Ref. 27) 129 437.304 1.4430 = 
Co-PhTT  [Co(C,H;N,S)2] 7.7 (Ref. 27) 103 413.37 1.7266 = 
Co-DMAET Na,[Co(C,H,oNS)(OH)2].2H,0 8.6 (Ref. 27) 103 283.366 1.2990 _ 
Co-TLA Na,[Co(C,H,0,S).2H,0] 8.9 (Ref. 27) 53 245.085 2.2503 — 
K,[Co&(CN) 0] 14.5 (Ref. 27) 77 533.73 5.7220 a 
[Co(NH 3)4(CO3)].NO3 10.8 (Ref. 27) 68 247.087 2.1184 a 
[Co(NH 3)4(NO,)2].NO3 11.7 (Ref. 27) 76 281.092 2.3364 co 
[Co(NH3)3(NO2)s3] 15.1 (Ref. 27) 64 248.06 3.2298 =? 
Cu-Complexes 
CuO 3.87+0.6 (Ref. 28) 3 79.54 6.8699 6.40 
CuSO,.5H,0 13.4 (Ref. 27) 71 249.69 2.481 2.284 
Cu-TDSA  H,[Cu(C,,Hs0,S).4H20] 8.0 (Ref. 27) 133 417.95 1.3376 = 
Cu-PPD Cu,[(C, ,H22N2S,)Cl4.4H20] 7.5 (Ref. 27) 154 587.317 ‘1.5381 = 
Cu-AcAcH Cu(C;H;0,), 12.3 (Ref. 27) 79 261.762 2.1298 ae 
Cu-AcAcCl Cu(C;H,0,C)), 10.9 (Ref. 27) 91 330.660 2.0610 = 
Cu-AcAcBr Cu(C,H,O,Br)2 10.4 (Ref. 27) 91 419.578 2.4960 = 
AcAcNO, Cu(C;H.O,NO,)2 8.8 (Ref.27) 106 351.762 = 1.5354 = 
Cu-TLA H,[Cu,(C,H,0,S)3].3H,O 15.6 (Ref. 27) 144 563.085 3.2833 = 
. the calculated and experimental values agreed fairly TDSA, Co-PhTT, Co-DMAET, Co-TLA, 
well. Thus Eq. (4) can also be applied for dielectrics. In K,[Co(CN):i0], [Co(NH3)4(CO3)JNO3, 
the present work, the calculated values of the plasmon [Co(NH3)4(NO2)2]NO3, [Co(NH3)3(NO2)3], Cu- 


energy for Co and Cu have been taken as 11.19 and 
10.79 eV respectively from our earlier papers!?"**. 


With the help of Eqs (1) and (4), we have obtained the 
following expression for the specific gravity: 
7 —— W-[5E+(ho,)y)? 
‘of ead 
Wolo pt 2829. 

; where suffix C stands for complexes or compounds. 
5 Eq. (6) help us in calculating the specific gravity of 
any complex or compound from its known chemical 
i - shifts SE. Using this equation, we have first calculated 
’ the specific gravities for the compounds CoS, CoS, 
 €oC,0,4, CoSO,4.7H20, CuO and CuSO,.5H,0 
whose specific gravities are known. Within the limits of 
_ the experimental error of chemical shifts, the calculated 
values for these compounds agree fairly well with their 
vn values of specific gravity. The present study 
tes that plasmons do exist in the above insulating 
and their energy can be calculated by the 
lence electron theory”®. This would mean 
a frequency is much larger than that 
the band gaps in these insulating 
fore, Eq. (6) has been used to calculate 
‘the complexes Co-PPD, Co- 


.. (6) 


TDSA, Cu-PPD, Cu-AcAcH, Cu-AcAcCl, Cu- 
AcAcBr, Cu-AcAcNO, and Cu-TLA whose specific 
gravities are unknown. 

In these calculations, the effective value of Z for 
cobalt and copper has been taken as one!*!*. The 
effective values of Z for oxygen, sulphur, chlorine, 
bromine, hydrogen, carbon and nitrogen have been 
taken as 6, 6, 7, 7, 1,4.and 3 respectively. These values of 
Z are their usual valencies. 

The calculated values of the specific gravities for 
cobalt and copper complexes are listed in Table 1 from 
which it can be seen that the error between the 
calculated and observed values for the specific gravity 
is in between 0.5 and 8%. This, as mentioned earlier, 
may be due to the experimental error which is of the 
order of 1eV in the measurement of chemical shifts. 

In conclusion, we would like to point out that the 
specific gravities of unknown compounds and 
complexes can be determined from their chemical 
shifts. This method, however indirect it may be, has the 
advantage of being free from the effects of macroscopic 
porosity and thus gives intrinsic specific gravity. The 
accuracy of the values of the specific gravity calculated 
by this method will depend upon how accurately the 
chemical shift is determined. 
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INSTRUCTIONS TO AUTHORS 


SCOPE 

The journal welcomes, for publication, full papers and 
short notes, reporting significant new results of research, in 
all areas of physics except space physics. The applied fields 
covered are electronics, electrical engineering, in- 
strumentation and applied mathematics. However, papers in 
applied mathematics with emphasis on only derivation and 
proofs and having no direct physical significance, will not be 
considered. Review articles are not published normally. 


SUBMISSION OF MANUSCRIPT 

Manuscripts for consideration should be submitted, in 
duplicate, to Editor, Indian Journal of Pure & Applied 
Physics, Publications & Information Directorate, Hillside 
Road, New Delhi 110012. They should neither have been 
already published nor be under consideration elsewhere. 

Manuscripts should be in English and typewritten on only 
one side of good quality paper, in double space, with 
adequate margin on all four sides. One original and one 
carbon or photo-copy, each complete in all respects 
including abstract, illustrations, appendixes, etc. are to be 
submitted. 


a PREPARATION OF MANUSCRIPT 

Authors may consult recent issues of the Journal to 
familiarize themselves with the general style and practices 

q adopted in regard to the various elements of a paper. 


General 

Manuscript should be presented in as concise a form as 
possible. Good attention should be given to spelling and 
a grammar. In giving names of chemical compounds and 
- structures, abbreviations of units of measurements, symbols 
cis "gaa notations, the style and practices recommended by the 
___—-: TUPAP and IUPAC, should be followed. — 

Bi _—s« Frequently repeating combinations of words, e.g. electric 
field gradient (EFG), junction field effect transistor (JFET), 
stimulated Raman emission (SRE), should be abbreviated 
subsequently, indicating the abbreviated form in paren- 

, _ thesis, as shown, at the place of their first occurrence. 

Pages should be numbered consecutively and arranged in 
the following order: Title, authors’ names with their 
utional affiliations and abstract, along with relevant 
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a “footnotes whenever necessary (on a separate sheet); 
introduction; experimental details/theory/method/analysis; 
ee ‘results; discussion; conclusion(s); acknowledgement; re- 
~ ferences and appendixes. Tables, captions for figures (with 
~ legends) and appendixes should be typed on separate sheets 


attached at the end of the manuscript. 
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e title should be neither too brief/general nor 

sarily long. It should reflect the content of the paper 
the maximum advantage in indexing. If a 
of a general series, a specific subtitle, 
icular aspect of the work covered in the 
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for the paper, the broad subject 
1ld be classified in the contents 
recent numbers of the journal for 
(indicated by an asterisk on the 
ress for correspondence, 
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highlighting the principal findings and conclusions. It should 
be in such a form that abstracting periodicals can use it 
without modification 


Introduction 

Long and elaborate introduction should be avoided. It 
should be brief and state the exact scope of the study in 
relation to the present status of knowledge in the field. 
Literature review should be limited strictly to what is 
necessary to indicate the essential background and the 
justification for undertaking the study. 


Materials, methods, apparatus, etc. 

The sources of materials and their purity, methods of 
preparation, procedure for measurements and their 
accuracies, etc. should be clearly stated to enable any other 
worker to repeat the work if necessary. New methods, 
techniques, theories, etc. should be described in adequate 
detail; but if they are well known, a mere literature reference 
to them will do; differences from standard ones, 
improvements or innovations should, however, be clearly 
mentioned. 


Results 

Only such primary data as are essential for understanding 
the discussion and main conclusions emerging from the 
study should be included. All secondary data as are of 
interest to a specific category of readership should not be 
included in the paper. Such data should be retained by the. 
authors for supply, on request, to any interested research 
worker. A footnote to this effect may be inserted at the 
relevant place in the paper. 

The results must be presented in a coherent sequence in a 
unified logical structure, avoiding repetition or confusion. 
Limitations of the results should be clearly stated. 

The same data should not be presented in both tabular and 
graphic forms. Only such tables and figures as are essential 
should be included. Simple linear plots that can easily be 
discussed in the text, should not be included. Infrared, 
ultraviolet, NMR and other spectra, DTA curves, etc. should 
be included only if they pertain to new compounds and/or 
are essential to the discussion; otherwise only significant 
numerical data should be given in the text or in a table. 


Discussion 
Long rambling discussion should be avoided. The 


discussion should deal with the interpretation of results 
without repeating information already presented under 
results. It should relate new findings to the known and 
include logical deductions. A separate section on 
‘conclusions’ can be given only when they are well 
established and of outstanding significance. Mere 
observation of qualitative trends of results should be 
distinguished from firm conclusions. Also, limitations, if 
any, to the conclusions should be clearly pointed out. _ 


Mathematical portions | 
Special attention should be given to the mathematical 


portions of the paper. Equations must be well separated from 
the text and written clearly with good separation between the 
successive lines. The usual norms of breaking long 
mathematical expressions should be adhered to. Equations 
should be numbered consecutively in Arabic numerals with 
the number in parenthesis near the right hand margin. 
‘Superscripts and subscripts should be clearly indicated in 
“per Bea and A sign respectively. Capital and small letters, 
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particularly of the same letter when both occur, as Well as 
letters or symbols likely to be confused one for the other, 
should be clearly distinguished. Special characters (e.g. Greek, 
script, vector, tensor, etc.) required must be indicated by 
marginal notes. Letters and symbols which should appear in 
bold face must be clearly indicated. To simplify typesetting: 
(i) long and complicated mathematical expressions which are 
frequently repeated should be replaced with single letter/ 
symbol, without clashing with the others used in the paper; 
(ii) the “exp” form of complex exponential functions should 
be used; and (iii) to simplify fractions, the solidus (/) is to be 
used and fractional exponents are to be used instead of root 
signs, e.g. 

write exp{ —i@(t, — t2)/2} and nat e~ ‘01 ~'2)/? 


2 j= we 2)1/2 40 K 5a 3, 
write (4@,; K3,/@ K p)’* and not ® K2 


Tables 

Tables should be numbered consecutively in Arabic 
numerals and should bear brief titles. Column headings should 
be brief. Units of measurement should be abbreviated and 
placed below the headings. Nil results should be indicated 
and distinguished clearly from absence of data. Inclusion of 
structural formulae inside the tables should be avoided as far 
as possible. Tables should be referred to in the text by 
numbers and not by terms like ‘above’, ‘below’, ‘preceding’ 
or ‘following’. Results should not be presented to a greater 
accuracy than that of the method employed. 


Illustrations 

The number of illustrations should be kept to the 
minimum. Wherever possible, e.g. a number of individual 
analogous figures referring to different variables, substances, 
molecules, etc. may be combined into one composite figure. 
All illustrations should be numbered consecutively in Arabic 
numerals. Captions and legends to the figures should be self- 
explanatory. Line drawings should be made with Indian ink 
on white drawing paper/cellophane sheet/tracing cloth, and 
drawn to approximately twice the printed size. 

The lettering should be uniform, preferably in stencil, so as 
to be not less than 1.5 mm after reduction widthwise to full 
page size (165 mm) or column size (80 mm). The size of 
geometrical shapes (used to distinguish different graphs), 
dots, lines, etc. should be sufficiently /arge to permit the 
necessary reduction without loss of detail. In the case of 
photographs, prints must be on glossy paper and contrast y. If 
an illustration is taken from another publication, reference 
to the source should be given and prior permission secured. 
Illustrations should be referred to in the text by numbers and 
not by terms like ‘above’, ‘following’ etc. 


Acknowledgements should not be exuberant and must be 
made only to real assistance rendered inconnection with the 
work reported in the paper. 


References 


should be numbered consecutively, in the order of their first 
occurrence, and should be indicated by superscript Arabic 


placement of references on numerals or other symbols should 
be avoided; in such cases the reference may be given in 
parenthesis in running text, e.g. “this yielded for na value of 
2.3 (Ref. 5)”. Full bibliographic details for all the references 
mentioned in the text should be listed in serial order at the 
end of the paper. 


In citing references to research papers, names and initials 
of authors should be followed, in order, by the title of the 
periodical in the abbreviated form (underlined), the volume 
number (two lines underneath), the year within circular 
brackets and the page number [e.g. Chandra B P & 
Shrivastava KK, J Phys & Chem Solids (GB), 39 (1978) 939). 
For names of periodicals, the abbreviations followed by the 
Physics Abstracts should be used. For periodicals not 
covered by Physics Abstracts, the title abbreviations should 
be according to the Bibliographic Guide for Editors and 
Authors, 1974, published by the American Chemical Society, 
Washington DC, USA; additionally the country from which 
the journal is published should be given in parenthesis 
immediately after the title abbreviation. If a paper has been 
accepted for publication, the names of the authors and the 
journal (with volume number and year, if known) should be 
given followed by the words “‘in press” [e.g. Wahi P K & 
Patel N D, Can J Spectrosc (Canada), in press.]. 
In references containing up to four authors, the names of all 
the authors with their respective initials should be given. The 
abbreviations et al., idem and ibid should be avoided. When 
there are more than four authors, only the names of the first 
three authors with their respective initials should be given, 
followed by the words ‘et al.’ 
Reference to a book should include details in the 
following order: name and initials of authors, the title of the 
book (underlined), name of publisher and place of 
publication within circular brackets and year and page (s) 
[e.g. Clayton G B, Operational amplifiers (Newnes- 
Butterworths, London), 4th Edn, 1977, 26]. If the reference is 
to the work of an author published in a book by a different 
person, the fact that it is cited from the source book should be 
clearly indicated [e.g Turnhout Van J, ‘Thermally 
stimulated discharge of electrets’ in Topics in applied physics: 
Vol. 33—Electrets, edited by C M Sessler Gpriaper Velae 
Berlin), 1980, 130]. 
Proceedings of conferences and symposia should be 
treated in the same manner as books. Reference to a paper 
presented at a conference, the proceedings of which are not 
published, should include, in the following order, names and 
initials of authors, title of the paper (underlined), name of the 
conference, and where and when it was held (e.g. Herczeg P, : 
Symmetry-violating kaon decays, paper presented to the ’ 
International Conference on High Energy Physics and 
Nuclear Structure, Vancouver, Canada, 13-17 August 1979). i 
Reference to a thesis should include the name of the 
author, title of the thesis (underlined), university or 
institution to which it was submitted and year of submission 
(e.g. Mehrotra S N, Many-bod; hnique. _ thei 
applications to interacting bosons, Ph D thesis, Ranch 
University, 1976). | aa 
Reference to a patent should include names of patentees, ae 
country of origin (underlined) and patent number, ‘the 
Organization to which the patent has been assigned (within | 
circular brackets), date of and 
reference to an a 
Labes M M, US Pat. 4,066,567 (to Temple niversity), 
January 1978; Chem. Abstr., 88 (No. 20) (1978), 138350 n]. 


j 


Dies: & Tans, alk ia 
Fibres & Pulps,_ Food: 


ak other Porest fa 0% ‘Name ‘in 
| languages, and { trade | names seg! provi 


t a L 


are pengienee: in detail Data pee area. oe 
and yield and import and export are provided. 4 


ee eae rene fications of raw materials, 
Be ee aS. well as finished products and statistical data 


‘ {7 le ous a apa q are ee The | regarding ‘production, demand, exports, 
:. articles are well illustrated. Adequate literature | imports. prices, etc., are provided: The artcle: 


| references are provided. . a ave suitably iipastea a Pi poacty- to the 
a Eleven volumes of the series covering letters | sources of information are pom: 
Fé A—Z have been published. | . 


Vol. WA- B) Rs. 80.00; Vol. IT (C) Rs. 95.00; Vol. 111 (D-E) | 
* Rs.,105.00; Vol.1V (F-G) Rs.65.00; Vol. IV: Suppl.|Fish & 
Fisheries Rs, 40.00; Vol. V(H-K) Rs. 75.00; Vol. VI(L-M) | i a; 
Rs. 90.00; Vol. VI: Suppl. Livestock Rs. 60.00; Vol. VII | Part 1(A-B) Rs. 54.00; Part II(C) Rs. 64.00; Part II (D-E) 

(N-Pc) Rs. 30.00; Vol. VIII (Ph-Re) Rs. 86.00; Vol. IX | Rs. 25.00; Part 1V(F-H) Rs. 25.00; Part V (I-L) Rs. 30.00, 
-(Rh-So) Rs. 104.00; Vol. X (Sp-W) Rs. 152.00; Vol. XI | Part V1(M-Pi) Rs. 28.00; Part VII (PI-Sh) Rs. 60.00; Part 
ee Rs. 102.00. VIII (Si-Ti) Rs. 66.00; Part IX (To-Z) Rs. 80.00. 


“Nine voltae’ of the series covering letiets a 
—Z have been Published. 


= 
ve 
“Ee te 


“HINDI EDITION: BHARAT KI SAMPADA—PRAKRITIK PADARTH ” ine sea 


WS : : a heen 
- x 


ia ies: oe Oe s ait to VI and two a ee of Wealth of nie Raw | hig eS ae oa, 
Materials series in Hindi already published. | Ps gee 
lish »d Volumes: ats Supplements: — tae ie = 


1 (33H) Rs. 38 Vol. I (#) Rs. 36: Vol. IIT (#7) Fish & Fisheries (Matsya & Matsyaki) Rs. a 
ae pone one aoe is Co as a Liverstock (Pashudhan aur Kukkut Palan) Rs. 34 


_ Vols. VII to XI under ane 


| Please contact: 
pet & Advertisement) 


PUBLICATION 
Proceedings:seminar.on primary —_ 
communications in Science 


7 


dophy C 


Birrese tt | 
Indian Thysanopter. 


30.00 13.00 
3.00 17.00 


emperatul 
a 


